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obvious  that  these  data  can  be  analyzed  and  interpreted  in  a  number  of  ways. 
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ABSTRACT 


^Spectral  and  cross- spectral  analyses  of  horizontal  wind,  components 
from  data  ~olle cted  by  the  author  on  an  experiment  at  the  upper  Texas  "oast 
are  presented.^.  ffl^as^rlrieErt  wa»~part*-of--ti»'-Se*r-areese  experiment  of  1968 
conducted, -by  the  Atmospheric  Science  Group  of  The  University  of  Texas  at 

*  V- 

Austin  .'-''The  wind  data  are  from  instruments  mounted  at  1,  6.7,  and  27  meter 
levels  on  two  towers.  One  tower  was  90  meters  from  the  coastline  and  the 
other  approximately  4.8  kilometers  inland  on  a  lin >  perpendicular  to  the 
coastline. 

The  data  were  recorded  simultaneously  at  both  towers  on  analog  tape 

and  digitized  electronically.  Data  samples  analyzed  are  from  a  continuous 

1 4- hour  experiment  which  covered  the  land  and  sea  breeze  regimes  of  the  sea 

breeze  phenomenon,  including  the  transition  period.  The  sampling  periods  are 

approximately  17,  67,  and  133  minutes,  (  '< 

An  autocorrelation  curve  and  Bpectrum  are  presented  for  each  data 

sample.  The  spectrum  is  analyzed  assuming  hypotheses  of  "white  noise"  and 

"red  noise."  Confidence  limits  are  placed  on  the  spectral  estimates.  The 

physical  interpretations  of  the  results  are  postulated  using  G.  I.  Taylor’s 

hypothesis.  Tor  the  nighttime  land  breeze,  it  was  found  that  the  spectra  of 

the  longitudinal  or  downstream  component  data  at  the  lower  levels  showed  more 

energy  at  the  higher  frequencies  than  that  of  the  lateral  component.  ThiB  was 

postulated  as  being  caused  by  mechanical  turbulence.  The  lateral  component 

indicated  long  meandering  waves.  Furthermore,  an  elongation  of  the  eddies 

with  the  wind  was  detected  using  the  67  minute  sampling  periods.  Ihe  size 

of  the  eddies  was  found  to  increase  upward.  The  cross-correlation  of  the 

Bame  components  between  levels  showed  higher  values  laterally  than  downstream. 

For  the  daytime  samples,  the  variance  increased  at  the  higher  frequencies  as 

iii 


solar  radiations!  heating  Increased.  The  downstream  elongation  was  not  detected 
as  local  horizontal  isotropy  was  approached;  furthermore,  the  horizontal  size 
of  the  eddies  did  not  increase  appreciably  with  height.  Th$  cross-correlation 
or  the  erne  ccspcssrts  between  levels  showed  high  correlation  values  laterally 
and  longitudinally.  From  the  data  and  the  cross -correlations,  it  was  indicated 
that  the  veering  of  the  winds  from  a  land  to  a  sea  breeze  regime  occurred  at 
the  top  level  first  and  at  the  bottom  level  last. 

The  spectra  of  the  daytime  onshore  flow  for  the  two  locations  are 
compared.  The  distribution  of  the  energy  among  the  frequency  ranges  is 
discussed  for  the  three  levels,  lie  loss  of  energy  by  the  longer  waves  and 
increase  in  energy  by  the  shorter  waves  as  convection  increases  are  discussed 
for  each  level  at  both  locations. 
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I.  INTRODUCTION 


S 

;■  A,  Statement  of  the  Problem 

T  .  The  geophysical  phenomenon  known  as  turbulence  has  been  the  subject 

1  of  numerous  investigations  which  have  led  to  theoretical  and  empirical  trea- 

tises  by  geophysicists,  mathematicians,  statisticians  and  others  in  related 
disciplines.  The  problem  of  explaining  the  nature  of  turbulence,  like  many 
others  in  physics,  can  be  approached  in  many  ways.  Scaoe  meteorologists  tend 
to  view  atmospheric  turbulence  as  a  microscale  phenomenon,  since  this  is  the 
most  noticeable  manifestation;  however,  departures  from  mean  flow  even  on  a 
global  scale  in  the  atmosphere  are  also  manifestations  of  turbulence.  In  the 
final  analysis,  these  and  intermediate  scales  of  turbulence  are  interrelated; 
therefore,  this  attempt  to  describe  the  aspects  of  a  particular  experiment 
of  wind  measurement  should  have  a  bearing  on  other  scales  of  turbulence. 

The  main  problem  dealt  with  in  this  particular  investigation  is 
that  of  determining  and  analyzing  the  spectra  of  the  horizontal  wind  velocity 
in  the  first  thirty  meters  of  the  atmospheric  boundary  layer. 

The  particular  hypothesis  to  be  examined  is  that  there  are  changes 
in  the  turbulent  structure  of  the  wind  caused  by  earth  surface  characteristics, 
thermal  stability,  and  radiational  cooling  and  heating;  and  furthermore,  that 
these  changes  can  be  detected  by  spectral  analysis.  The  description  of  the 
physical  phenomenon  then  follows  from  the  interpretation  of  the  spectra. 

The  wind  data  are  from  wind  instruments  mounted  at  three  discrete 
elevations  on  two  separate  towers.  The  towers  vere  situated  approximately  4.8 
kilometers  (3  miles)  apart  on  a  line  perpendicular  to  a  coastline t  one  tower 
about  90  meters  inland  from  the  surf  and  the  other  4.8  kilometers  (km)  inland. 

The  location  was  on  the  Texas  costs t  approximately  65  km  northeast 
of  Galveston,  Texas  -  an  area  characterized  by  a  long  straight  coastline  and 
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flat  terrain  consisting  of  marshland  and  grass .  This  van  an  ideal  location 
for  observing  the  wind  from  the  Gulf  as  it  is  modified  by  the  land  mass. 
Different  types  of  flov  were  also  obtained  because  of  the  land-sea  breeze 
phenomenon. 

The  experiment  van  part  of  a  major  meteorological  undertaking  by 
the  Atmospheric  Science  Group  of  the  University  of  Texas  at  Austin  during  the 
months  of  June  and  July,  1968.  The  experiment  is  described  in  Chapter  EL. 
Various  phases  such  as  the  instrumentation,  calibrations  procedures,  and  data 
reduction  are  covered  in  detail. 

The  analyses  and  discussion  of  the  data  sure  given  in  Chapter  III. 
Data  samples  include  offshore  flow  or  land  breeze,  and  onshore  flow  car  sea 
breeze,  with  other  data  samples  covering  the  transition  periods.  Prior 
studies  of  the  sea  breeze  phenomenon  on  a  mesoscale  have  been  made  in  the 
same  area  by  Hsu  (1967,  1969)  and  by  McPherson  (1968).  The  study  presentd 
here  is  a  microscale  analysis. 

The  spectral  and  cross-analyses  and  the  interpretations  are  also 
given.  The  representations  are  discussed  in  detail.  The  results  of  the 
analysis  are  summarized  in  Chapter  XV. 

Prom  the  viewpoint  of  data  gathering,  this  experiment  utilized 
s ante  of  the  latest  techniques  in  the  recording,  digitizing,  and  data  re¬ 
duction  phases.  The  method  of  observation  is  one  of  the  foundations  upon 
which  the  final  analysis  depends j  in  this  experiment,  the  use  of  electronic 
tape  recorders  reduced  the  observer  error  and  eliminated  the  inherent  time- 
lag  response  error  of  a  mechanical  chart  recorder. 

One  year  after  this  particular  experiment,  another  project  of  a 
much  larger  scale  was  undertaken  by  several  federal  governmental  agencies, 
universities,  private  industry,  and  National  Center  for  Atmospheric  Research 
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in  the  Barbados  Island  region,  with  the  cooperation  of  the  Barbados  govern¬ 
ment  ( Kuettner  and  Holland,  1969)*  Hie  data- gathering  techniques  included 
some  similar  to  the  one  described  in  this  experiment.  With  the  new  Global 
Atmospheric  Research  Program  (GAHP)  in  the  planning  stages  (Chamey,  1969), 
undoubtedly  other  methods  will  be  found  to  measure  the  physical  phenomena  of 
the  atmosphere  and  to  interpret  the  results. 

B.  Atmospheric  Turbulence 

1.  Statistical  Description 

The  geophysical  problem  of  turbulence  In  the  atmospheric  boundary 
layer  may  be  most  conveniently  approached  by  considering  the  earth's  surface 
as  a  rigid  boundary  which  is  heated  during  the  daytime  be  solar  radiation  and 
is  cooled  at  night  by  infrared  radlational  heat  loss.  This  diurnal  cyclic 
process  in  turn  heats  and  cools  the  adjacent  air  which  flows  over  the  earth's 
surface.  The  immediate  effect  is  the  density  gradient  and  associated  positive 
or  negative  houyancy  created  by  the  temperature  differences  resulting  in  the 
rising  or  lowering  of  the  air.  At  this  point  other  factors  such  as  wind 
velocity,  roughness  of  terrain,  and  gravity  become  important. 

In  the  laboratory,  Osborne  Reynold’s  familiar  classical  experiment 
of  the  stream  of  dye  flowing  through  a  pipe  showed  that  laminar  flow  existed  . 
until  a  certain  value  of  the  Reynolds  number  was  achieved.  Reynolds  number 
is  defined  by  Re  =  Ud/v,  where  U  is  the  mean  velocity  of  the  fluid,  d  is  the 
diameter  of  the  pipe,  and  v  is  the  kinematic  viscosity,  i.e.,  the  ratio  of  the 
dynamic  viscosity  ’  ~>  the  density.  Above  this  critical  value,  the  character  of 
the  flow  becomes  "disorderly"  or  "turbulent",  causing  the  velocity  profile, 
which  was  a  parabolic  conic  section  during  the  laminar  state,  to  change  its 
shape.  The  single  streak  of  the  dye  is  completely  distorted  downstream.  This 
phenomenon,  which  may  be  called  mixing,  cai  _.es  an  ever-increasing  transfer  of 


heat  and  momentum  downstream  from  one  layer  to  adjoining  layers.  According  to 
Reynold s ,  the  composite  nature  of  the  flow  Is  one  of  a  mean  motion.  This  com¬ 
plicated  motion  is  of  a  random  character;  a  specific  "eddy" or  fluctuation 
cannot  he  predicted  in  space  or  time  by  an  explicit  mathematical  formulation. 

The  concept  of  a  mean  motion  with  an  eddying  motion  superimposed  on 
it  may  perhaps  he  restated  as  a  stochastic  geohpyslcal  process .  This  may  he 
interpreted  as  a  geophysical  phenomenon  made  up  of  a  deterministic  component 
and  a  non-deterministic  component,  the  mean  motion  being  the  deterministic 
part.  A  deterministic  process  is  defined  as  one  which  may  be  expressed  by  an 
explicit  mathematical  function,  e.g.,  a  sine  wave.  A  completely  non-determi- 
nistic  process  represents  complete  chaos,  e.g.,  a  random  process.  In  general, 
all  geophysical  processes  are  stochastic  in  nature,  in  the  sense  that  even 
orderly  processes  usually  contain  a  non-determlnistic  component. 

Applying  this  concept  to  the  atmosphere,  a  someviiat  laminar  type  of 
flow  may  be  envisioned  during  the  nighttime  hours  under  clear  skies  and  stable 
thermal  statification  in  layers  near  the  ground;  however,  because  of  the  surface 
roughness  and  variations  in  the  density  of  the  medium,  the  layers  immediately 
adjacent  to  the  ground  (lower  boundary)  are  not  laminar.  This  differential 
density  process  is  not  present  in  the  (Reynolds)  laboratory  experiment  of  the 
dye  streak  and  Increasing  fluid  velocities.  The  problem  then,  albeit  analogous, 
is  more  complicated  in  the  atmosphere  even  under  uiie  most  favorable  atmospheric 
conditions . 

Notwithstanding,  the  problem  may  be  attacked  from  a  physical  point 
of  view  and  the  equations  of  motion  of  an  incompressible  viscous  fluid  and  the 
equation  of  continuity  may  be  used  as  a  foundation. 

If  a  geophysical  realization  such  as  an  anemometer  recording  of  the 
wind  velocity  component,  u,  in  the  x  direction  of  a  Cartesian  coordinate  system 
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is  given,  the  realization  may  be  said  to  represent  a  stochastic  process.  Thus 
u  may  be  represented  by  a  mean  motion  plus  a  random  or  fluctuating  motion,  or 
u  «  u  +  u1 .  The  over bar  represents  a  mean  value  defined  by 

T/2 

u  «  lim  5  f  u(t )  dt 
T  -*  w  T  -T/2 

where  T  Is  the  time  interval  of  the  realization  and  t  is  the  time.  The 
fluctuation  Is  represented  by  u1 .  (in  order  for  the  u’  to  be  random,  the 
mean  value  of  u*  must  be  equal  to  zero.)  Similarly,  the  component  in  the  y 
and  z  directions,  v  and  w,  respectively,  are  given  by  v  *  v  +  v*  and  w  * 
w  +  w* . 

The  turbulent  form  of  the  Navier-Stokes  equations  of  motion  may  be 
written  in  an  Eulerian  frame  as 

(t  -  p  u2  -  p  u'2)  + (t  -  p  u  v  -  p  u’v*) 

xx  K  y  '  by  xy 

(\x  "  P  “  ^  -  P  «T*7) 

'v-’"-'™  *4  ( v ■ p yS ■ p 'r,2) 

(t  -  p  v  w  -  p  v’v*  )  (l«l) 

-  p  U  w  -  p  u'wf)  +  -g|  (T  -  P  V  w  -  p  v’v* ) 

~2  ~2, 

(rzz  -  p  w  -  p  v*  ) 


o  a  A  p  +  ~ 

K  ot  K  x  dx 


dz 


P  ST  =  p  Fy  +  to 


dz 


dw 

dt 


=  p  P  +  r2 
v  z  dx 


dz 


5 


*  £  StVs-'iswr  *gr  • 


where  F  ,  F  ,  F  are  the  external  forces  such  as  friction,  gravity,  Coriolis, 
x  y  z 

per  unit  mass,  a  is  density  of  the  medium,  and  r  ,  r  ,  t  and  t  ,  t  ,  t 

'  w  J  xy*  xz’  ya  xx’  yy  za 

are  the  mean  tangential  and  normal  viscous  stresses.  The  Reynolds  stresses 

""""2  _ _ _  1 . . -  2 

(also  called  apparent  or  virtual  stresses)  are  p  u'  ,  p  u'v*,  p  u'v' ,  p  v*  , 
p  w'2,  and  p  v'w’.  These  express  the  effect  of  the  velocity  fluctuations  in 
transporting  momentum  across  a  specified  surface  in  the  fluid. 

It  is  important  to  note  that  these  stresses  are  in  the  f^rm  of 
correlation  coefficients,  for 


JLZ-1L. 
Ns  s 
x  y 


are  the  equations  for  the  correlation  between  variables  x  and  y,  where  8^  and 
s  are  the  sample  standard  deviations.  This  shows  that  the  virtual  stresses 

y 

are  proportional  to  the  correlation  coefficient  between  the  components  of 
turbulent  velocity  observed  at  a  fixed  point  in  an  Eulerian  frame  of  reference. 

G.  I.  Taylor's  theorem  (l92l)  showed  that  the  dispersion  of  a  group 
of  particles  in  a  statistically  uniform  field  of  flowj  i.e.,  one  in  which  the 
mean  eddy  energy  p  u’  is  constant  with  respect  to  time  and  space,  is  deter- 

'  2 

mined  by  the  magnitude  of  the  velocity  variance  of  the  fluctuations  u*  and 
the  nature  of  the  correlation  between  fluctuations  which  act  on  the  particle 
at  different  times. 

Specifically,  the  Lagrangian  correlation  is  given  by 


R(S)  =  u'(t)  u^t  + 
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and  denotes  the  correlation  of  u*  at  times  t  and  t  +  g. 
Integrating  Eq.  (1.2)  with  respect  to  £, 


u*2  J  R(  5)  dg  =  J  u’(t)  u'(t  +  |)  a§ 


u'(t)  J  u’(t  +  5}  df 


u’(t)  X 


(1.3) 


where  X  is  the  distance  traveled  by  the  particle  in  a  time  t  because  of  the 
eddy  motion.  Since 


u‘"(t)X  ~  X2  (lA) 

the  following  equation  may  be  written! 

_ -  _ _  t 

X2  =  2  u*2  J  J  H(§)  d§  dt  (1.5) 

o  o 

This  shows  that  the  distance  travelled  In  time  T  by  the  particle  is  expressed 

p 

by  T,  u1  and  R(§).  Hence  the  diffusion  of  the  particle  is  reduced  to  a  sin¬ 
gle  quantity  —  the  correlation  coefficient  between  the  particle  velocity  at 
time  t  and  time  t  +  5. 

If  the  time  of  diffusion  t  is  small  so  that  the  correlation  is 
essentially  still  close  to  the  initial  value  of  unity,  then  by  Eqe.  (lA)  and 

(1-3) 


1 
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and 


(1.6) 


I f  the  diffusion  is  in  a  stream  of  air  moving  with  velocity  U  and  the 
dispersion  is  observed  at  a  distance  x  downstream  from  the  source,  then,  by 
substituting  t  =  x/U  in  Eq,  (1.6) 

7  x2  7  u*2 

<  '■  ■  ■  ■■  =X  »'■■■  — . . 

x  U 


The  analysis  also  yields  a  length  L^,  given  by 


X  u‘  -  1^.  7  u‘2  =  u'2  J*  °  R(5)  d? 


or 


=  7  u»2  r* 


J  R(l)  <15 


(1.7) 


when  the  irregular  motion  is  such  that  a  time  t  is  found  in  which  R(5)  =0 
for  all  §  greater  than  tQ.  This  length,  L^,  is  called  the  Lagrangian  inctegral 
scale  (length).  According  to  Taylor  (1935,  P>  425),  this  "...bears  the  same 
relationship  to  diffusion  by  turbulent  motion  that  the  mean  free  path  does  to 
molecular  diffusion."  Taylor  (op.  cit.)  showed  that  the  theory  of  the  "scales" 
of  turbulence  could  be  applied  to  either  the  Lagrangian  or  the  Eulerian  con¬ 
cepts  of  fluid  flow.  His  method  was  first  used  by  Mr.  L.  P.  G.  Simmons 
(Taylor,  1935,  p.  444j  1937,  p.  499)  who  used  a  wind  tunnel  to  find  experi¬ 
mentally  the  correlation  between  the  turbulent  components  of  velocity  u  and 
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u  at  two  points  a  distance  y  apart  across  the  direction  of  flow  behind  a 

y 

honeycomb  with  square  wire  mesh  grid.  Statistically,  if  the  correlation  R(y) 
is  plotted  against  y,  the  curve  should  represent  the  distribution  of  u  along 
the  y  axis.  If  the  points  are  many  "eddy”  diameters  apart,  the  correlation 
should  drop  to  zero.  At  this  point  of  zero  correlation,  an  integral  "scale" 

y 

is  defined  by  Lg  =  J  R(y)  dy,  referred  to  by  Taylor  {1935,  p.  ^26)  as  the 
o 

"average  size  of  the  eddies"  in  an  Eulerian  system. 

Thus,  the  two  scales  defined  are 


y 

^  =  J  R(y)  dy  (length) 

o 

I u'2  J*  °  R(§)  d§  (time) 
u  o 

where  ^  is  the  Lagrangian  integral  scale  of  Eq.  (1.7). 

Later,  Simmons  and  Slater  (lay lor,  1938)  showed  how  the  energy  of  the 
fluctuations  could  be  distributed  among  the  motions  of  different  scale  lengths 
or  eddy  diameters,  i.e.,  the  existence  of  a  spectrum  of  turbulence .  In  his 
article  Taylor  (op.  cit.)  discussed  the  "connexion  between  the  spectrum  of 
turbulence  measured  at  a  fixed  point,  and  the  correlation  between  simultaneous 
values  of  velocity  measured  at  two  points."  This  is  explained  in  the  following 
section. 

2.  Spectrum  of  Turbulence 

The  plots  of  x  and  y  wind  components  shown  in  Fig.  1.1  are  examples 
of  actual  values  of  the  raw  data  obtained  in  the  field  experiment  to  be  des¬ 
cribed  in  Chapter  II.  Note  that  the  fluctuations  appear  to  be  randomj  i.e., 
there  are  no  discemable  periodic  components  in  the  data. 
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Pig.  1.1  Plots  of  the  x~  and  y-components  of  Velocity  at  1 
Meter  Level  on  the  "Beach"  lover,  1047-1104  CUT 


If  the  x-component  of  eddy  velocities  is  resolved  at  a  fixed  point 

g 

into  harmonic  components,  the  variance  u*  is  the  sum  of  the  contributions 

from  all  frequencies.  This  value  u'  is  a  measure  of  the  average  level  of 

turbid  ent  energy  over  all  frequencies  involved  in  the  turbulent  motion. 

Between  the  frequencies  f  and  f  +  df ,  there  is  some  energy  fraction. 

"I"1'  V)  Vj 

Therefore,  if  the  total  energy  is  given  by  u’  ,  than  u*  P(f)df  is  the  energy 
2 

contribution  to  u’  from  frequencies  between  f  and  f  +  fd.  Thus 

J*  F(f)df  =  1. 

J  o 

Also, 


u,2=/  u’2  P(f)df  (1.8) 

o 

A  form  of  the  spectrum  curve  is  obtained  if  the  spectrum  function  F(f)  is 
plotted  against  frequency  f. 

By  definition,  Fourier  integrals  are  definite  integrals  which  repre¬ 
sent  an  arbitrary  function  in  an  unlimited  interval  (Carslav,  1950>  P-  311 ) . 

O 

To  show  that  u*  is  the  sum  of  all  contributions,  the  Fourier  integral  and  a 
form  of  Parse  veil's  theorem  are  used. 

Let  u*  =  0  (t).  The  Fourier  coefficients  are  given  by 

i,  (f)  -  A 
<r>  ■  5 


J  0  (t)  COS  2  TT  ftdt 

03 

J  0  (t)  sin  2  tt  ftdt 


(1-9) 
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and  therefore, 


j  t  —  ^ 


(t)  *  2  TT  OOS  2  TT  ft  +  I2(f)  SllS  2  TT  ftj  df  (l.lO) 


Parseval’s  Theorem  states s 


I ,r  >2) 


Applying  the  theorem  to  Eq.  (1.10), 

J  [♦  (t)]2  dt  =•  2  n2  J*  [^  +  ^]  df 


(1.11) 


(1.12) 


If  the  integrals  in  Eq.  set  (1.9)  and  the  left  side  of  Eq.  (1.12)  are  taken 
over  a  long  period  of  time  T  instead  of  infinity,  then,  for  a  constant  variance 
of  velocity  fluctuations, 


♦  (tr 


2  M 


/  liffi 
'  o  T  -* 


4  *4 


df 


(1.13) 


Comparing  this  result  with  the  turbulent  spectrum  of  Eq.  (l.  8),  ve  can  write 

u'2  F(f )  =  2  tt2  lim 

T  • 

2  4  +  4  -2 

where  the  quantity  2  tt  lim  — - -  ^  ■  —  is  the  contributions  to  u*  from  the 
components  of  frequency  between  f  and  f  +  df. 

At  this  point,  a  general  interpretation  of  the  spectrum  curve  and 
the  correlation  curve  obtained  from  a  geophysical  realization  such  as  the  afore¬ 
mentioned  (Fig.  1.3)  plots  of  the  x-  and  y-components  of  wind  velocity  will  be 
given.  Two  general  cases  are  considered}  namely,  (l)  variations  from  small 
eddies  carried  by  the  wind  stream  (of  mean  velocity  U)  past  a  fixed  point,  and 
(2)  variations  from  large  eddies  carried  by  the  wind  stream. 


4  +  4 


d.i^) 
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For  the  first  case  of  small  eddies,  the  energy  should  be  confined 
more  to  the  high  frequencies;  therefore,  the  spectrum  curve  should  show  higher 
energy  contributions  In  the  higher  frequencies  than  In  the  case  of  the  larger 
eddies.  However,  the  correlation  function  should  be  just  the  opposite  in  that 
the  case  of  the  longer  eddies  show  higher  correlations  at  longer  lag  (time) 
periods  than  in  the  case  of  the  small  eddies.  This  is  easy  to  conceive,  since 
the  fluctuations  in  the  this  case  are  more  rapid.  Also,  this  should  result  in 
high  correlation  a  short  distance  x  from  the  fixed  point  and  decreasing  corre¬ 
lation  as  x  increases  away  from  the  point  of  measurement. 

Taylor  (1938,  p.  478)  assumed  that  if  the  wind  velocity  which  carried 
the  eddies  was  very  much  greater  than  the  turbulent  velocity,  the  changes  in 
u‘  at  the  fixed  p-.  int  were  due  to  an  unchanging  pattern  of  turbulent  motion 
over  the  point,  or 


u*  =  u’(t‘ )  =  u«U/5) 

where  U  is  the  mean  wind  velocity.  Note  that  the  ratio  u’/u  ->  0  as  U  increases . 
Assuming  this  is  true  for  small  values  of  the  ratio,  then 


R(t)  =  v.fo) 

.2 


and 


R(t)  =  u'.fo)  u'  (*■  .+ 

„»2 
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Generalizing  Barseval’s  Theorem  and  applying  it  as  in  Eq.  (1.13)*, 


u*(t)  u'(t  +  x/0)=  2tt2  f  11m 

o  T  •* 


X2  +  I2 

1  a 


COB  { 2  TT  P  —  )  dLf 


and  from  Eq.  (l.lk). 


V 


(1.15) 


u'(t)  u*(t  +  x/0)  *»  u’2  f  P(f)  COB  (2  tt  F  —  ) 

"o  U 


(1.16) 


Thus 


alLtl-a^JUs/g).  mf  r (f)  CQ8  (2  w  Pi  )  at 

,,.2  °  0 


(1.1T) 


or 


R(x)  «  f  F(f)  cos  (2  tr  F  —  )  df 

‘  o  ir 


(1.18) 


By  Fourier  transforms  the  Fourier  integral,  for  f(x)  an  even 
function,  may  he  written 


f(x)  »  p2-.  f  g  (p)  cos  pxdp 
V  2  n 


g(p) 


2 

J  f(x)  coe  pxdx 
V  2  IT  O 


(1.19) 


♦This  discussion  follows  essentially  Haurwitz  (1968), 


1U 


Substituting  the  following  values  in  the  above  equations  i 


p  «  f  F(x)  a  R(  x) ,  and  g(p)  * 


U  F(f) 
2  V2tT 


Eq.  set  (1.19)  becomes 


R(x)  -  J  F(f)  cos  (2  tf  F  ~  )  df 


U 


F(f)  -if  R(x)  cos  (2  tt  F  -  )df 
U  o  U 


(1.20) 


Therefore,  R(x)  and  the  spectrum  function  ■■••p£=2  are  Fourier  transforms  of 

2v2  nr 

one  another]  that  is,  if  the  correlation  function  of  the  values  at  time  u  and 
at  time  t  +  t'  is  known,  the  spectrum  function  may  be  found,  and  vice-versa. 

3.  Scales  of  Turbulence 

The  two  scales  derived  by  Taylor  were  the  Iagrangian  integral  scale 

given  byL.  =*  u*  J*  R(§)dg  and  the  Eulerian  scale  given  byL£  -  J  R(y)dy. 

o  o 

The  Lagrangian  scale  was  derived  as  an  anology  to  the  mean  free  path  of  mole¬ 
cular  motion;  the  Eulerian  length  scale  was  interpreted  as  a  possible  "average 

size  of  the  eddies."  Another  scale  which  may  be  obtained  is  the  Lagrangian 

00 

time  integral  scale  given  by  t^  *  J*0  R(5)^5. 

The  correlation  functions  computed  from  the  data  obtained  throughout 
the  experiment  to  be  described  in  Chapter  III  were  such  that  the  Eulerian  time 
integral  scale  is  given  by 


*E  =  Jo  R(t)dt 


(1.21) 
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The  existence  of  such  an  integral  scale  Is  sufficient  to  imply  an  ergodlc 
process .  This  Is  based  on  the  equality  of  the  ensemble  average  and  the  tempo¬ 
ral  average  of  any  member  of  the  ensemble  -which  may  be  found  as  follows i 
Given  a  Bt&tionary  ensemble  of  stationary  functions  u(t),  the  ensemble  average 
of  u(t)  is  the  expected  value  S  j  u(t)  |.  The  time  or  temporal  average  la 

T/2 

lim  £  f  u(t)  dt 
T  ■+  to  '-T/2 


The  condition  that  these  averages  be  the  same  may  be  met  by  requiring  that 
the  ensemble  average  of  the  square  of  the  difference  between  the  two  averages 
converge  to  zero  j  that  is. 


-1C*  C  -  *  I  u(t)  I  j2 1  -*  0  as  T  ->  » 


(1.22) 


Lumley  (1962)  has  shown  that  a  sufficient  condition  for  the  expected  value  to 
converge  to  zero  is  that  the  Integral  obtained  by  expanding,  interchanging 
integration  and  ensemble  averages,  and  noting  that  the  correlations  R(t)  »  R(-t) 
from  Eq.  (1.22) 


f  (1  -  |)  R(t)  dt  -*  0  as  T  -»  •  (1.23) 

’  o 

where  R(t)  is  the  correlation  coefficient,  converge  to  a  finite  value.  If 

eo 

J*^  R(t)  dt  »  tg  exists,  which  is  some  finite  time  interval,  then  the  integral 
In  Eq.  (1.23)  has  a  finite  value. 

The  Fourier  transform  pair  of  the  Eulerian  space  correlation  R(x) 
and  the  spectrum  function  F(f)  given  in  Eq.  set  (l.2o)  were 
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R(x)  -  JQ  F(f)  cob  ZJLJLx  df 

F(t)  r  B(x) 

U  o  U 

From.  Taylor's  hypothesis,  t  »*  x/5,  "where  t  is  the  time,  x  is  the  distance 
between  two  points  and  U  is  the  mean  wind  speed 


R(t)  -  f  F(t)  cos  2  n  Ftdt 
’  o 

09 

Ht)  -  4  J*  R(t)  cos  a  Tt  PMt 
o 


(I.!*) 


00 

At  frequency  zero,  F(o)  *  h  f  R(t)dt,  and  from  Eq.  (l.2l),  F(0J  =  4t£.  mils 

'  o 

states  that  the  Eulerian  time  integral  scale  is  equal  to  one-fourth  the  zero 
intercept  of  the  spectral  estimate.  This  zero  frequency  corresponds  to  a  wave 
of  infinite  length  (a  straight  line),  which  implies  again  that  the  integral 
scale  must  exist  and  be  finite)  otherwise  the  eddy  energy  in  the  atmosphere 
would  be  infinite  (Munn,  1966  p.  76). 

Panofaky  et  aL  (1958)  have  shown  that  R(t)  and  R(x)  are  nearly  equal 
for  low  levels  of  turbulence  intensity,  which  is  given  by  the  ratio  of  the 
standard  deviation  of  the  wind  speed  and  the  mean  downstream  wind  velocity. 

The  relationships  between  the  Iagrangian  and  Eulerian  special  and  temporal 
correlations  are  more  complex,  and  only  a  few  studies  of  these  have  been  made. 
Corrsin  (1963)  found  that  the  Iagrangi*»r  time  integral  scale  is  roughly  equal 
to  the  Eulerian  length  integral  scale  divided  by  the  root-mean-square  velocity 
at  high  Reynolds  numbers.  Similarly,  Eulerian  time  scales  give  Values  approxi¬ 
mately  equal  to  the  Iagrangian  scales. 
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A  basic  method  of  viewing  the  relationships  is  to  follow  Taylor's 
hypothesis.  The  lagrangian  time  scale  t^  is  given  by  ^/u  where  L ^  is  the 

lagrangian  length  scale  and  u  is  the  eddy  velocity}  the  Eulerian  time  scale  is 
given  by  tg  * L  g/u  where  1  g  is  the  Eulerian  length  scale  and  U  is  the  mean 
velocity.  The  ratio  of  to  should  be  inversely  proportional  to  the  in¬ 

tensity  of  the  turbulence,  I.  This  may  be  written  as 


U  I  t„ 


E 


Pas quill  (1963)  also  showed  that  the  Eulerian  autocorrelation  R(tg)  is  equal 
to  the  Lagranglan  autocorrelation  R(t^)  when  g  ■  6  tg. 

The  basic  theory  upon  which  the  spectrum  analysis  of  the  wind  compo¬ 
nents  is  based  has  been  presented  in  the  above  discussions.  Interest  now  can 
be  centered  on  two  main  areas}  namely,  (l)  obtaining  the  data  to  be  analyzed, 
and  (2)  applying  the  necessary  statistical  operations  and  recent  developments 
in  spectrum  analysis  to  get  meaningful  information. 
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II.  EXPERIMENTAL  PROCEDURES 

A,  Location  and  Description  of  Site 

The  experiments  were  conducted  by  the  Atmospheric  Science  Group  of 
the  University  of  Texas  (Austin)  on  the  'Texas  Gulf  coast  approximately  forty 
miles  northeast  of  Galveston,  Texas  during  the  month  of  June,  1968.  Pig.  2.1 
is  a  map  of  the  area.  The  coastline  is  nearly  straight  for  more  than  forty 
miles j  inland,  the  terrain  is  flat,  rising  only  approximately  ten  feet  in 
fifteen  miles  perpendicular  to  the  coast. 

Two  telescoping  towers  were  erected,  each  30  meters  in  height  when 
fully  extended  (manufactured  by  the  Rohn  Co.  y  j  one  was  located  approximately 
90  meters  from  the  beach  and  the  other  4.8  km  inland  on  a  line  perpendicular 
to  the  coast  near  the  intracoastal  canal  Bhown  in  Fig.  2.1.  (Henceforth,  the 
towers  will  be  referred  to  as  the  "Beach"  tower  and  the  "Canal"  tower.)  The 
only  vegetation  between  the  two  towers  consisted  of  different  types  of  grasses 
about  one-third  to  nearly  one  meter  in  height.  The  Canal  tower  was  erected  on 
a  clearing  about  10  meters  in  diameter.  Otherwise,  there  were  no  structures 
or  obstacles  to  interfere  with  the  wind  flow  between  and  around  the  two  towers. 

B.  Description  of  Instruments 

1.  Wind  sensors!  The  wind  direction  and  i/ind  speed  sensing  instru¬ 
ments  utilized  thoughout  the  experiment  were,  respectively,  Gill  microvanes 
and  Gill  3-cup  anemometers  manufactured  by  the  P..  M.  Young  Co.  Some  of  the 
specifications  (nominal  values)  furnished  by  t-he  manufacturer  are  as  follows j 

Microvane:  Model  12301  -  delay  distance  (50$  recovery),  0.95  ® 
damping  ratio,  0.44 

maximum  overshoot  for  sinusoidal  fluctuations,  25$ 
construction  -  aluminum  alloy  with  foa**  polystyrene  fins 
dimension  -  72,5  cm 

3-cup  anemometer)  Model  12101  -  distance  constant  (63$  recovery)  =  2.8 
threshhold  sensitivity  with  tachometer  generator  0.45  - 
0.54  mps 
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The  power  supply  translators  were  built  by  Mr.  Tom  Phillips  of  the  Atmospheric 
Science  Group  of  the  University  of  Texas. 

The  time  constant  of  a  geophysical  transducer  may  be  defined  as  the 
period  required  for  the  transducer  to  respond  to  63.2  percent,  £(l  -  l/e)> 
where  e  ks  2 . 718  J  of  a  step  change  in  the  signal.  In  this  case,  the  geophysical 
phenomenon  is  wind;  therefore,  a  more  appropriate  term  "distance  constant"  is 
used  to  describe  the  nominal  performance  of  the  instrument.  The  distance  con¬ 
stant  of  a  wind  sensor  is  defined  as  the  length  of  air  column  required  to 
cause  the  sensor  to  respond  to  63.2  percent  of  the  square-wave  change  in 
speed.  This  term  more  aptly  describes  the  performance  because  it  is  indepen¬ 
dent  of  wind  velocity;  e.g.,  if  the  time  constant  were  3.0  seconds  for  an 
instrument  taken  from  still  air  to  an  environment  of  3.0  mps,  then  the  same 
instrument  would  have  a  time  constant  of  only  1.5  seconds  if  It  were  taken 
from  still  air  to  a  6  mps  wind  speed.  In  each  case  the  same  length  of  air 
must  pass  for  the  stepwise  change  of  63*2  percent.  J'rom  the  specification 
given  above  for  the  anemometer,  a  wind  speed  of  2.71  mps  would  give  63.2 
percent  change  in  one  second. 

2.  Temperature  Sensors 

The  wet  bulb  and  dry  bulb  temperatures  were  measured  by  bead 
thermistors  housed  in  a  device  developed  by  the  Atmospheric  Science  Group 
at  the  University  of  Texas.  The  wet  bulb  temperature  sensor  was  covered  with 
an  elongated  muslin  wick  which  had  one  end  immerBed  in  distilled  water. 

Q 

The  thermistors  were  accurate  to  +  0.1  Celsius  and  had  a  time  con¬ 
stant  of  approximately  one  second.  A  squirrel-cage  fan  was  used  to  continuously 
aspirate  the  system  internally  where  the  sensors  were  mounted. 

A  complete  sounding  was  taken  every  eight  minutes.  This  slow  ascent 
or  descent  of  the  device  gave  ample  time  for  the  sensors  to  adjust  to  the 


environment  and  record  representative  readings  since  the  ascent-descent 
velccity  was  about  5  cm  per  second. 

3.  Recorders 

The  D.  C.  signals  of  the  wind  direction  and  wind  speed  sensors  were 
recorded  on  eight  channel  Precision  Instrument  analog  tape  recorders,  Model  6100, 
One  recorder  was  used  at  each  of  the  sites.  The  number  of  available  analog 
recording  channels  was  the  limiting  factor  on  the  number  of  wind  measurement 
levels  recorded.  Simultaneous  recordings  of  wind  speed  aid  direction  were 
obtained  for  the  1  m,  6.7  m,  and  27  m  levels  at  both  towers.  One  of  the 
channels  was  used  to  record  reference  pulBes,  the  other  for  voice  recordings 
of  time  and  tape  footage  for  reference.  All  recordings  were  made  at  0.375 
inches  per  second,  allowing  nearly  l4  hours  of  continuous  analog  recording  on 
a  standard  1800-foot  reel  of  magnetic  tape. 

The  wet  bulb  and  dry  bulb  temperatures  were  recorded  on  Esterline- 
Angus  Graphic  Ammeter  strip  chart  recorders.  Reference  marks  for  time  and 
tape  footage  coincided  with  the  voice  recordings  on  magnetic  tape.  These  strip 
chart  recorders  also  show  the  position  of  the  instrument  package  by  a  separate 
side  marking  pen  for  the  top  and  bottom  points  of  the  temperature  package 
cycle . 

The  tape  recorders  and  strip  chart  recorders  wore  installed  in  air- 
conditioned  van- type  vehicles.  The  air  conditioning  and  fans  were  installed 
to  cool  and  circulate  the  air  inside  the  vehicle  as  a  precaution  against 
malfunctions  of  the  tape  recorders.  This  was  a  necessary  measure  since  the 
site  was  located  in  humid  marshland  and  the  experiments  were  conducted  in 
June. 

C.  Instrumented  Towers 

The  wind  directional  vanes  and  3“CUP  anemometers  described  in  Section 
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B  were  mounted  on  the  Beach  tower  on  2  meter  long  aluminum  booms  at  1*0,  1*9, 
3*0,  4*5,  6.7,  11*7,  and  27.0  meters  above  the  ground.  At  the  Canal  tower  the 
instruments  were  mounted  at  1.0,  3.0,  6.7,  11.7.  and  27.0  meters  on  5  seperate 
booms.  The  booms  were  placed  on  the  towers  so  that  any  interference  induced 
by  the  tower  structural  members  would  coincide  with  the  approximate  5-degree 
electrically  open  zone  of  the  potentiometers  in  the  wind  vane  as  shown  in 
Jig.  2.2.  Any  expected  interference  would  probably  occur  when  the  wind  was 
from  a  north-northeasterly  direction. 

Only  three  levels  were  used  for  the  turbulence  measurements  described 
in  this  report;  these  were  the  1.0,  6.7,  and  27,0  meter  levels. 

The  dry  and  wet  bulb  temperature  sensors  described  in  Section  B 
were  mounted  on  a  rope-pulley  mechanism  which  raised  and  lowered  the  sensors, 
completing  one  cycle  every  16  minutes  from  1.5  meters  to  25-3  meters. 

D.  Instrument  Calibration 

The  Gill  Microvanes  were  calibrated  in  the  laboratory  at  the 
University  of  Texas  at  Austin  prior  to  the  experiments.  Calibration  curves 
for  microvane  #1  is  shown  in  Fig.  2.3.  The  calibrations  were  repeated  Just 
prior  to  mounting  the  instruments  to  insure  that  no  damage  had  occurred  in 
the  packing,  transporting,  and  unpacking.  The  wind  speed  calibrations  of  the 
3-cup  anemometers  made  in  the  laboratory  and  prior  to  conducting  the  exper¬ 
iments  were  limited  to  a  single  maximum  output  calibration  at  1800  BFM  (equal 
to  22,35  mps)  listed  in  the  specification.  The  maximum  output  voltages  for 
tape  recordings  were  set  at  1.000  volt  (+  0.005),  which  equaled  ^  9  mps,  since 
no  speeds  higher  than  this  value  were  anticipated.  After  performing  the  ex¬ 
periments,  the  anemometers  were  checked  again  for  the  calibration  value  at 
1800  RIW,  Wind  tunnel  calibration  measurements  at  the  low  wind  speed  values 
were  made  using  a  hot-wire  anemometer  as  a  cross-reference.  The  regression 
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Fig. 


2.2  Top  View  of  Lowest  Section  of  Telescoping  "Beach" 
Tower  Showing  Orientation  of  Tower  and  Location 
of  Aluminum  Boom 


Calibration  Curve  for  a  Microvane 


curves  were  computed  for  each  of  the  microvares.  and  each  of  the  anemometers 
using  a  least-squares  technique.  Table  II. 1  shows  the  regression  equations 
for  the  six  sensors  used  in  the  tower  experiment.  The  equations  are  of  the 
form  y  =  mx  +  b,  where  b  is  the  y  intercept,  m  is  the  slope,  and  x  the  abscissa. 

All  data  pertaining  to  a  particular  sensor  were  reduced  using  its 
respective  calibration  curve. 


Table  II.l—Regression 
Equations  for  Wind  Sensors 

Micro  vane  No.  3-  cup  Anemometer 

(x  in  degrees  of  arc,  (x  in  meterr  per  second, 

y  in  volts)  y  in  volts) 

y  =  0.1555x  -  0.01115 

y  =  0.1587x  -  0.02388 

y  =>  0.1577x  -  0.01625 

y  =  0.l632x  -  0.03967 

y  =  0.1516x  -  0.01067 

y  =  0.1607x  -  0.02661 


1.  y  =  0.004008x  -  0.1050 

2.  y  =  0.004074x  -  0.1005 

3.  y  =  0.004028x  -  0.1171 

4.  y  =  0.003911x  -  0.0843 

5.  y  *  O.OO3976X  -  O.0975 

6.  y  =  0.004026x  -  0.0995 


E.  Analog-to-Digltal  Conversion 

The  time  series  of  signals  recorded  in  the  field  experiments  were 
digitized  in  the  laboratory  using  the  equipment  shown  in  the  block  diagram, 

Fig.  2.4. 

The  two  other  speeds  were  available  for  playback,  namely  3  3/4  and 
7  l/2  ips.  The  analogue  to  digital  voltmeter  (Dana  model  5403)  coupled  with 
the  digital  tape  recorder  (Precisian  Instruments  Model  II67)  had  a  sampling 
rate  which  could  be  varied  up  to  40  samples  per  second.  A  low  frequency 
oscillator  and  trigger  generator  (square  wave)  were  used  to  control  the  sample 
rate;  these  are  shown  as  the  sample  rate  generator  in  the  block  diagram. 

In  order  to  reduce  the  time  of  digitizing,  a  playback  speed  of  3  3/4 
ips  was  used  with  a  sampling  rate  of  5  samples  per  second.  Since  this  speed 
is  10  times  faster  than  the  recording  speed,  the  real  time  was  100  seconds  of 
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Pig.  2,5  Example  of  Aliasing 


digitized  data  for  10  seconds  of  digitizing  time.  At  a  rate  of  5  samples  per 
second,  the  real  time  sampling  rate  equaled  0.5  samples  per  second,  or  2  seconds 
per  sample.  This  is  the  sampling  used  throughout  the  entire  set  of  digitized 
data  samples,  although  other  sampling  rates  were  tried  to  test  for  aliasing. 

The  speed  of  2  seconds  per  sample  was  found  to  he  the  optimum.  The  Speed omax 
W  recorder  (Leeds  and  Northrup  Co.)  shown  in  the  diagram  as  the  Chart  Recorder, 
served  "both  as  a  secondary  signal  monitor  and  as  a  quick  data  record.  With  it, 
the  performance  of  the  tape  recorder  could  be  monitored,  and  the  geophysical 
signal  plot  served  as  a  reference  for  actual  visual  interpretations  of  the 
sigials  and  permanent  records.  Interpretations  of  stationarity,  best  intervals, 
and  instrument  behavior  were  valuable  as  first  approximations  by  use  of  these 
charts. 

The  problem  of  aliasing  is  always  present  when  analog  records  are 
digitized.  The  digitized  data  or  data  sample  are  said  to  be  aliased  If  under¬ 
sampling  occurred  in  the  digitizing  of  the  data.  Fig.  2.5  shows  an  example 
of  aliasing.  A  signal  with  frequency  of  4  Hz  is  shown  by  the  solid  line.  The 
circles  depict  the  sampling  points,  giving  an  apparent  frequency  of  1  Hz  as 
shown  by  the  dotted  line. 

Mathematically,  if  h  is  the  sampling  Interval,  the  Nyquist  frequency 
is  given  byf N  =  l/2  h.  This  frequency  is  the  highest  frequency  that  may  be 
detected  with  data  sampled  at  intervals  h.  All  other  frequencies  higher  than 

f  will  be  "folded"  into  a  lower  frequency  between  zero  and  l/2  h. 

N 

There  are  essentially  two  ways  of  avoiding  this  aliasing  problem. 

One  is  to  determine  the  sampling  interval  h  prior  to  digitizing)  the  other  is 
to  pre- filter  the  raw  data  to  remove  the  frequencies  higher  thanfR.  The  first 
method  has  been  used  throughout  all  the  samples  uBed  in  the  experiments. 

Another  point  to  be  considered  is  that  the  time  series  obtained 
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ver»  from  physical  instruments  which  had  inherent  inertia.  These  time  series 
are  said  to  be  "Band- limited”  for  no  frequencies  higher  then  a  certain  maximum 
frequency  f^  are  contained.  This  frequency  is  determined  by  the  frequency 

response  of  the  instrument.  This  feature  than  facilitates  the  task  of  finding 
the  Nyquist  frequency  because,  for  a  given  set  of  atmospheric  conditions  to 
be  measured,  “  fjj»  This  frequency  may  be  found  (approximately)  by  using 
the  specifications  and  calibration  data  of  the  instrument. 

Another  problem  to  be  considered  initially  is  one  of  whether  to 
make  analog  recordings  or  to  digitize  directly  from  the  sensor,  for  the  latter 
is  less  time  consuming  and  more  direct.  The  problem  involves  the  proper 
sampling  rate;  undersampling  causes  aliasing  problems  and  oversampling  is 
inefficient.  The  matter  of  convenience  of  the  direct  digitizing  method 
sacrifices  the  optional  sampling  rate  which  is  the  main  advantage  of  the 
analog  method.  Per  example,  10  to  l4  hours  of  recordings  might  be  made  of  a 
geophysical  process  such  as  the  wind.  Afterward,  the  digitizing  interval 
might  be  varied  to  test  for  aliasing  e’-d  to  find  the  optimum  rate.  If  con¬ 
siderable  changes  occur,  such  as  wind  speed  of  2  mps  gradually  increasing  to 
10  mps,  the  optimum  sampling  rate  will  differ  in  each  case,  and  can  be  readily 
changed.  On  the  other  hand,  for  a  direct  digitized  record,  the  sampling 
rate  or  rates  have  been  fixed,  and  cannot  be  varied.  Constant  monltering 
by  the  investigator  is  needed  when  the  data  are  being  recorded  in  order  to 
change  the  sampling  rate  if  required. 

P.  Data  Reduction 

The  analog  records  on  magnetic  tape  were  first  played  back  »,tvI 
recorded  on  strip  charts  using  the  Speedomax  recorder.  Time  "hacks"  and 
footage  were  noted  on  the  strip  chartB  for  reference  marks .  One  set  of  data 
from  &  tower  thus  required  six  different  strip  chart  recordings.  Appendix  A, 
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Fig,  A.1,  shows  a  sample  frees  the  beach  location  for  June  11-13,  1968.  The 
records  were  than  examined  for  any  instrument  malfunctions  and  power  failures 
which  may  have  occurred.  More  important,  however,  was  the  actual  first 
approximation  of  any  obvious  trends.  Also,  since  wind  statistics  were  to  be 
computed,  different  sample  domains  could  be  obtained  by  visual  observation 
of  the  record.  For  example,  the  gross  statistics  such  as  second  moments  about 
the  (variance)  would  be  different  for  stable  atmospheric  conditions  as 

compared  with  unstable  conditions ,  for  two  different  populations  are  involved . 
Although  the  statistics  could  be  computed,  bimodal  distribution  and  other 
configurations  could  result.  By  a  visual  inspection  of  the  data,  "best" 
sampling  periods  may  be  obtained. 

Once  the  sampling  periods  of  records  were  determined,  the  analog 
records  were  digitized  by  playing  one  channel  at  a  time  through  the  DANA 
digital  voltmeter.  The  sampling  interval  was  determined  by  means  of  an 
electronic  clock  mechanism  as  ehown  in  Fig.  2.4.  The  reference  pulses  on 
one  of  the  channels  served  as  starting  "marks"  far  the  digitizing  of  all  the 
channels,  since  the  starting  points  had  to  be  exact.  Strip  chart  recordings 
were  aleo  made  during  the  digitizing  process.  An  oscilloscope  was  also  used 
as  a  visual  monitor  of  the  output  signal,  since  the  Sraedooax  recorder  had 
an  inherent  time  delay  and  high  frequency  interference  could  not  be  detected 
as  readily  as  on  the  oscilloscope. 

The  number  of  data  points  in  each  sampling  record  was  between  4200 
and  6400j  thus  a  continuous  record  without  inter-record  gaps  or  end-of-files 
was  obtained  for  8400  seconds  (l40  minutes)  up  to  12,800  seconds  (213  minutes). 

The  half-inch  magnetic  tapes  with  the  discrete  data  samples  were 
then  processed  in  the  Control  Data  Corporation  66 00  computer  at  the  University 
of  Texas  Computation  Center-  Additional  checks  were  made  on  the  quality  of 
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the  data. 


These  were! 


(1)  The  digital  voltmeter  occasionally  quantized  a  value  of  9-999 
volts.  In  a  record  of  500  samples,  about  4  of  these  9*999 
values  would  be  found.  Rather  than  compute  a  statistical 
analysis  of  the  adjoining  points  to  substitute  for  these 
high  values,  the  two  adjoining  values  were  simply  added  and 
divided  by  two.  The  result  was  substituted  for  the  high  value 
in  question. 

(2)  The  number  of  data  points  selected  for  each  data  sample  was 
4000.  The  data  sample  of  4000  was  taken  from  the  4200  to  6400 
data  points  at  a  '’best"  interval  to  insure  some  degree  of 
stationarity  for  the  data  samples. 

(3)  Certain  segements  of  the  data  were  plotted  and  compared  with 
the  strip  chart  records  to  check  for  aliasing.  The  raw  data 
were  then  converted  into  wind  direction  in  degrees  and  wind 
speeds  in  meters  per  second.  The  individual  calibration  curves 
obtained  by  the  linear  regression  best  fit  were  used  for  each 
of  the  data  points  from  a  particular  sample.  The  vaxaes  ob¬ 
tained  for  wind  speed  and  direction  were  printed  out  for 
further  visual  inspection. 


III.  DATA  ANALYSIS  AND  DISCUSSION 
A.  Data  Samples 

The  data  samples  for  the  study  are  taken  from  a  continuous  experi¬ 
ment  which  began  at  2200  CUT  on  the  12th  of  June  and  ended  on  the  13th  at 
1300  CDT,  This  continuous  run  helped  insure  that  the  instruments  and  equip¬ 
ment  all  were  calibrated  the  same  since  daytime  and  nighttime  comparisons  -were 
to  be  made.  The  data  samples  and  the  code  to  be  used  throughout  this  chapter 
are  as  shown  in  Table  IH.1  for  evening  samples  and  Table  III. 2  for  the  day¬ 
time  samples. 

The  evening  samples  all  occurred  during  a  land  breeze.  The  differences 
between  them  are  the  sampling  periods  and  the  time  of  observation.  The  daytime 
samples  are  more  diversified  in  that  the  sea  breeze  phenomenon  exists  during 
the  perl  .d  encompassing  the  samples,  first  at  the  beach  and  later  inland. 

The  land-sea  breeze  circulation  ideally  occurs  under  certain  condi¬ 
tions,  one  of  which  is  that  during  this  period  the  area  surface  winds  should 
not  be  j "trolled  by  a  synoptic  flow  pattern.  During  the  evening  when  the  land 
surface  cools  by  radiation,  the  sea  surface,  because  of  Its  physical  charac¬ 
teristics,  maintains  a  more  constant  temperature.  Therefore,  the  air  will  flow 
from  land  to  sea.  This  land  breeze  is  characterized  by  very  light  winds, 
usually  less  than  5  mps. 

During  the  morning  daylight  hours,  as  the  land  surface  becomes 
warmer  than  the  sea  surface,  the  flow  is  reversed.  This  reversal  of  the  wind 
flow  is  usually  characterized  by  veering  (clockwise  turning)  of  the  wind  with 
time.  The  daytime  samples  were  taken  during  the  veering  of  this  wind  and  30 
minutes  after  the  wind  veered  from  ENE  to  ESE.  Hsu  (1969)  made  a  mesoscale 
study  of  the  sea  breeze  using  data  obtained  previously  with  surface  observations, 
radiosondes,  pilot  balloons,  and  aircraft  observations .  A  study  of  the  land 
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Table  III.l.  Evening  Cases  Investigated 
All  land  Breeze  Occurrences 


Code 


Beach  Night 


Canal  Night 


Date  and 
Time  (car) 


Period  Sampled 
( sec) _ 


I-BN 

I-CN 

12/2213-2230 

1000 

II- BN 

II- CN 

12/2213-2320 

4000 

III-EN 

III-CN 

12/2320-13/0027 

4000 

IV- BN 

IV- CN 

12/2213-13/0027 

8000 

Table 

III. 2.  Day  Cases 

Investigated 

Code 

Date  and 
Time  ( CDT) 

Period  Sampled 
(sec) 

Remarks 

(Beach  Day) 
I- HD 

13/1000-1017 

1000 

Easterly  Flow 

II-BD 

-3/IOI7-IO34 

1000 

Beginning  of  Sea  Breeze 

III-BD 

13/1005-1022 

1000 

Beginning  of  Sea  Breeze 

IV- BD 

13/1047-1104 

1000 

30  min  after  Sea  Breeze  Started 

V-BD 

13/1103-1120 

1000 

Same  time  as  Canal  IV- CD 

VI- BD 

13/1000-1107 

4000 

Beginning  of  Sea  Breeze 

VII-BD 

13/1103-1214 

4000 

Sea  Breeze  Regime 

(Canal  Day) 
I- CD 

13/1103-1120 

1000 

Land  Breeze  Regime 

II- CD 

13/1155-1212 

1000 

Beginning  of  Sea  Breeze 

III- CD 

13/1236-1253 

1000 

30  min  after  Sea  Breeze  Started 

IV- CD 

13/110^-1210 

4000 

Land  Breeze  Generally  (Same  time 
as  Beach) 

V-CD 

13/1155-1302 

4000 

Beginning  of  Sea  Breeza 

VI- CD 

13/1210-1317 

4ooo 

Sea  Breeze  Regime 

VII- CD 

13/1103-1317 

8000 

Land  and  Sea  Breeze  Regimes 

breeze  in  this  same  area  has  been  made  by  Beit  (1969),  McPherson  (1968)  des¬ 
cribes  the  sea  breeze  with  a  3-dimensional  numerical  model. 

In  order  to  perform  the  analysis  of  the  horizontal  wind  components 
the  digitized  data  in  the  form  of  wind  speed,  V,  and  wind  direction,  0,  were 
transposed  to  horizontal  components  X  and  Y  by  the  orthogonal  trans format ion 

X  *  J V J  sin  (6  +  tt) 

Y  =  |  Vj  cos  (6  +  it) 

The  X  is  positive  for  a  wind  from  the  west,  and  negative  for  an  east  wind}  Y 
is  positive  for  a  wind  from  the  south  and  negative  for  a  north  wind.  Thus, 
conventional  vector  notation  agrees  with  meteorological  convention  since  a 
resultant  vector  is  180  degrees  from  a  resultant  wind  vector.  The  statistics 
for  a  bivariate  normal  distribution  (Crutcher,  1957)  were  computed  for  these 
components;  a  graphical  representation  of  these  statistics  is  shown  in  Fig.  3.1 
for  a  hypothetical  case.  The  conic  section  with  a  vector  radius  equal  to  the 
standard  deviation  theoretically  contains  approximately  63  percent  of  the  (X,  Y) 
coordinates  or,  adding  180  degrees  to  the  vector  angles,  thiE  ellipse  contains 
63  percent  of  the  wind  origins.  The  ellipse  1b  obtained  from  the  various 
statistical  parameters  computed  as  shown  in  Appendix  D  for  the  x-  and  y-com- 
ponents . 

The  sample  records  then  were  grouped  into  500  data  points  each  by 
taking  500  samples  for  the  1000- sec  periods  (l  sample  per  2  seconds),  arith¬ 
metically  averaging  4  consecutive  samples  to  get  500  samples  for  the  4000  sec 
periods,  and  averaging  every  8  for  the  8000  sec  periods.  A  few  examples  of 
these  500-point  records  for  the  x-  and  y-components  are  shown  in  Appendix  B 


for  daytime  and  evening  cases.  Statistical  frequency  distributions  were  made 
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‘stWJtU.-A  T, 


Standard  Deviation  and  Angle  of  Rotation  of 
Major  Axis 


of  the  sample  records  and  a  chi-square  goodness- of-f It  test  for  normality  also 
was  made  (Appendix  C).  The  bivariate  statistics  computed  are  shown  in  Appendix 
D. 


Before  computing  the  spectra,  the  coordinate  axes  were  rotated  to 
coincide  with  the  resultant  or  mean  wind  vector  by  the  tranformatlons 

u  =  X  cos  6  -  Y  sin  S 
v  =  X  sin  9  +  Y  cos  9 


where  9  is  the  angle  of  rotation.  This  was  done  in  order  to  obtain  a  down¬ 
stream  wind;  that  is,  a  mean  wind  parallel  to  the  x-axis .  Taylor's  hypothesis 
could  then  be  applied.  Therefore,  the  downstream,  or  longitudinal,  wind  is 
the  u-component  and  Its  mean  value  is  the  mean  wind  speed.  The  cross  wind,  or 
lateral  wind,  is  the  v- component  and  has  a  aero  mean  value. 

B.  Spectrum  and  Cross- Spectrum  Analysis  Formulae 

In  Chapter  I,  the  development  of  the  statistical  theory  of  turbulence 
showed  the  relationship  between  the  autocorrelation  function  and  the  modified 
spectrum  function.  The  spectrum  analysis  technique,  as  performed  on  the  data, 
is  a  fairly  new  development.  According  to  Tukey  (1968),  the  beginnings  of 
modern  spectrum  analysis  were  developed  independently  by  him  and  by  Bartlett 
in  the  late  19^0* s.  More  recent  developments  include  the  utilization  of  the 
"fast  Fourier  transform"  (Gentlemen  and  Sande,  1966). 

The  following  formulae  are  essentially  those  of  Blackman  and  Tukey 
(1958)  and  Bendat  and  Fiersol  (1968).  Some  of  the  problems  involved  in  obtain¬ 
ing  the  spectral  lines,  such  as  aliasing,  stability,  truncation,  etc.,  have 
been  investigated  by  Eddy,  Duchou,  and  Almaz an  (1968). 


Autocovariance  Functions 


rrm  1  .  a  J-  -»  _ _  V.  ~  A 

i  t  t  1 B  lUlUJCAUli  wuc 
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value  on  the  neighboring  values.  Por  a  purely  random  series  the  values  are 
independent  and  no  correlation  is  found  j  however,  when  the  measurements  are  of 
a  geophysical  signal,  there  will  normally  exist  some  correlation  between  the 
values  of  x(t)  at  time  t  and  time  t  +  T,  where  r  is  the  lag.  The  autocovariance 
function  may  be  expressed  by 


C(r)  *  /  ,f  x(t)  x(t  +  t)  dx(t)  dx(t  +  t)  (3.l) 

-»  -co 

for  a  stationary  time  series.  Prom  the  ergodic  hypothesis,  ensemble  averages 
are  equal  to  time  averages,  hence  Eq,  (3.1)  may  be  written 

T/2 

C(t)  -  lim  s  /  x(t)  x(t  +  t )  dt  (3.2) 

T  -»  «  1  IT/2 

The  autocovariance  C(r)  and  autocorrelation  function  R(t)  are  related  by  the 
ratio  R(t)  =*  C(r)/C(o)  where  C(o)  is  the  covariance  at  lag  zero.  This  form  shows 
that  division  by  C(o)  is  equivalent  to  normalizing  the  covariance  because  the 
covariance  at  lag  zero  is  equal  to  the  variance.  All  the  raw  data  were 
initially  normalized  by  the  procedure  z  =  (x  -  x)/s  where  s  is  the  sample 
standard  deviation. 

The  form  of  the  summation  equation  used  in  the  computations  is 


R(t)  “  FTT  L.  *N  +  k  k  =  0,  1,  2,  ...,  M.  (3.3) 


where  M  is  the  maximum  lag.  The  value  for  the  maximum  lag  was  taken  as  10  per¬ 
cent  of  the  total  number  of  data  points. 

Spectral  Function  and  Estimate  1  The  raw  spectral  estimate  for  normal - 
lized  discrete  data  with  a  mean  of  zero  may  be  given  by 
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M-l 

Gx(f)  -  2  fa  [  R(0)  +  2  £  Rr  cos  (£-£-£)  +  Ry  eoe  (-5^)  ]  (3.4) 

where  h  is  the  sampling  interval,  fjj  the  Hyauist  frequency,  M  the  maximum  lag, 
and  r  a  o,  1,  2,  ...,  M.  The  raw  spectral  estimates  were  computed  for  the  nor¬ 
malized  data  by  using  the  following  form 

M-l 

Gx(f)  -  2  h  [  R(0)  +  2  £  Rr  cos  £-jp  +  (-l)k  Rjj  ]  (3.5) 

r=l 

where  k  =  0,  1,  2,  ...,  M,  the  harmonic  corresponding  to  the  frequency  f  =  kf^/M. 
This  gives  the  spectral  estimate  at  discrete  frequencies  for  the  positive  and 
negative  sides  of  the  two-sided  spectrum. 

Smoothing  of  Spectral  Estimates:  Ifcere  are  several  "windows"  which 
are  used  in  smoothing  spectral  estimates.  One  of  the  simplest  is  "Hanning", 
named  after  the  Austrian  meteorologist  von  Hann  who  smoothed  meteorological 
data  (hut  not  for  spectrum  analysis)  by  using  the  coefficients.  (Other  windows 
sometimes  used  by  spectrum  analyzers  are  Hamming,  Parzen,  and  Bartlett,  plus 
others.)  The  Hanning  lag  window  weighting  function  may  be  defined  in  terms 
of  a  multiplier  given  by 


D(k)  *  (l  +  cos  ~Jj“) 
=  0 


k  =  0,  1,  2,  . ..,  M 
k  >  M  (3.6) 


Using  the  spectral  estimates  of  the  data  obtained  by  Eq,  (3.5),  the 
equivalent  smoothing  of  the  estimates  wwre  computed  by 
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G*  -  0.5  G  +  0.25  0, 
o  o  l 

GX  “  0,5  Go  +  °*5  GX  +  0,25  G2 

Gk  "  °*25  Gk-1  +  °*5  Gk  +  0.25  Gk+1  k  «  2,  3, 

<£  -  °-5  W  +  °*5  °m 


(3.7) 


M 


Cross-covariance  Function  and  Estimate:  The  croBs-covar lance  func¬ 
tion  for  two  sets  of  data  x  and  y  is  a  measure  of  the  dependence  of  the  values 
of  x  on  those  of  y.  The  cross-covariance  estimates  for  a  data  sample  may  be 
taken  exactly  as  for  the  autocovariance  function  given  by  Eq.  (3.l)j  that  is, 


T/2 


C  (t)  lim  |  f  x(t)  y(t  +  t)  dt 

rj\  _  x  \_m/n 


-T/2 


(3.8) 


Similarly,  the  cross-correlation  function  for  the  sampled  data  was  computed 
using 

N-k 


and. 


Rxy  ^  a  N-k  Z  *N  yN  + 
N-l 


N-k 

v(T)  =A  I 

N-l 


(3.9) 


N-k  L  yN  *N  +  k 


Croos-Spectral  Function  and  Estimate!  The  cross-spectral  function 
for  two  time  series  is  the  Fourier  transform  of  the  cross-correlation  function, 
exactly  as  for  the  single  record  Fourier  transform  pair.  Tale  cross-spectral 
function,  however,  is  mathematically  a  two-sided  complex  valued  function 
given  by 


:xy^f'  =  LCVf'  "  1  Hxy1^ 


U.10; 
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where  Co^  Is  the  ”co~spectral”  function  (even)  and  iQ^,  the  complex  term,  is 
the  "quadrature"  function.  The  co-spectrum  is  a  measure  of  the  covariance 
between  the  in-phase  components  of  the  two  series}  whereas  the  quadrature 
spectrum  is  a  measure  of  the  out-cr- phase  frequency  components,  the  out-of¬ 
phase  being  one- quarter  period. 

Eq.  (3.10)  may  also  be  written 

-i«  (f) 

Gxy(f)  “  |Gxy<f)l  e  ^  <3.ll) 

where  the  modulus,  or  magnitude,  is  defined  as 

lGxy  (f)i  “{[Coxy{f)]2  +  [V(f)]2}l/2  (3*12) 


and  the  phase  angle  $  by 


<f>  ■  t“1  cVf)  /0v<f>] 


(3.13) 


The  magnitude  spectrum  shows  whether  the  "amplitude"  of  the  component  at  a 
frequency  f  in  series  x  is  associated  with  a  large  or  small  amplitude  at  the 
same  frequency  in  series  y.  The  phase  angle  determines  whether  the  frequency 
components  In  series  x  lead  or  lag  the  components  of  series  y  for  a  given 
frequency. 

In  the  data  analysis,  the  raw  cross -spectral  estimates  were  computed 
for  M  +  1  discrete  frequencies  at  k  harmonics  corresponding  to  f  =  fjjk/M  for 
k  =  0,  1,  2,  ...,  M  by 


Co. 


?  h 


rA 

L  O 


M-l 


9  ) 
—  / , 


cos 


-l)k  A..  ] 


r-1 


(3-14) 


4o 


where 


and 


where 


Finally, 


Ar  -  0.5  [Rxy(T)  +  ^(t)]  *  -  0,  1,  M 


M-l 


r»l 


y  “  0,  1,  2,  .  .«,  M. 


the  raw  estimates  were  Hanned  by  the  following: 


Co  (o)  =>  0.5  Co  (1)  +  0.25  Co  (2) 

Cd  (l)  -  0.5  Co  (1)  +  0.5  Co  (2)  +  0.25  Co  (3) 

Cd  (k)  =  0.25  Co  (k  -  l)  +  0.5  Co  (k)  +  0.25  Co  (k  +  1) 

Cd  (M)  =>  0.5  Co  (M)  +  0.5  Co  (M  -  l) 

Q*  (0)  »  0.0 

T 

Q  (l)  =»  0.5  Q  (l)  +  0.5  Q  (2)  +  0.25  Q  (3) 

t 

Q  (k)  =  0.25  (k  -  1)  +  0.5  Q  (k)  +  0.25  Q  (k  +  l) 

t 

Q  (M)  =  0.25  Q  (M  -  1) 


Coherence  Function:  The  coherence  function  may  be  given  by 


,  ,  lGxv  Wl 

Coh  (f)  = 


4l 


(3.15) 


(3.18) 


where  G  and  G  are  the  variance  spectral  estimates.  This  form  of  the  coher- 

*  y 

ence  function  Is  real- valued  and  has  a  magnitude  range  from  zero  to  one.  It 
is  a  measure  of  the  correlation  between  time  series  x  and  time  series  y  at 
particular  frequencies.  If  the  value  is  zero,  the  series  are  said  to  be  in¬ 
coherent  at  that  f requency _s  If  completely  coherent,  the  value  of  coherency  is 
unity. 

The  coherence  function  was  computed  using  the  Harmed  variance  spec¬ 
tral  estimates  for  digitized  data  given  by  Eq.  (3.7)  and  Hanned  co-spectral 
and  quadrature  estimates  given  by  Eqs.  (3.16)  and  ( 3 . 17 ) •  Thus, 

[Co'(f)]2  +  u'(f)]2 

Coh  (f)  =  — - - -  (3.19) 

G'  (?)  G'  (f) 

a  y 

Is  the  equation  used. 

Rejection  Filtration t  A  filter  usually  refers  to  an  algebraic  ex¬ 
pression  used  on  a  set  of  data  to  isolate  periodic  components  of  a  particular 
size  either  to  examine,  modify,  or  remove  these  components  from  the  set  (Eddy 
et  jlL,  Chap.  9)*  The  process  of  filtering,  when  used  in  conjunction  with 
variance  spectrum  analysis,  is  sometimes  referred  to  as  "prewhitening".  This 
implies  that  the  spectral  extimates  are  filtered  in  an  attempt  to  flatten  or 
make  more  even  the  spectral  estimates  so  as  to  resemble  those  of  the  "white 
noise"  spectrum  to  be  discussed  in  Sec.  C  below. 

As  in  nearly  all  aspects  of  spectrum  analyses,  indiscriminate  (not 
in  a  statistical  sense)  use  of  filters  may  result  in  erroneous  significant 
peaks  in  the  spectral  lines. 

The  need  for  two  spectral  analyses,  one  using  raw  data  and  another 
using  filtered  data,  is  evident  in  working  with  geophysical  data  involving 

+  «.««« „  « ,  +  r-\  TT  — - . .  1  ^ - 
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involved  the  mean  wind  and  fluctuations  superimposed  on  it.  This  mean  wind 
constitutes  the  longer  waves.  Intuitively,  if  the  raw  data  are  used  to  com¬ 
pute  the  variance  spectrum  for  a  stable  atmosphere  at  a  height  above  5  m, 
nearly  all  of  the  variance  will  be  concentrated  toward  low  frequencies.  The 
small  amount  of  variance  which  is  distributed  among  the  higher  frequencies 
will  be  difficult  to  discern.  Once  the  data  have  been  analyzed  to  assess 
the  importance  of  the  large  low  frequency  spectral  peaks,  the  data  may  be 
filtered  to  reject  the  low  frequencies.  This  allow  the  higher  frequencies  to 
be  examined  for  fluctuations  by  performing  a  second  spectrum  analysis  on  the 
filtered  data. 

In  the  actual  analysis,  spectral  estimates  were  obtained  using  the 
raw  data  (normalized)  and  the  results  examined.  Since  the  object  was  to  re¬ 
move  the  low  frequencies  for  the  second  spectral  analysis,  a  difference  filter 
given  by  y(t)  =»  x(t)  -  x(t  -  l)  was  used  on  the  original  wind  components. 

Hie  gain  function  of  this  filter  is  given  by 

Gain  (f)  -  2  jsin  rr  f| 

The  computation  and  graphical  plot  of  this  gain  function  is  shown  in  Appendix 
D.  This  spectral  window  multiplies  the  Hanned  spectral  estimates  by  (2-2cos2rtf) 
as  Bhovn  in  the  appendix;  therefore,  if  the  "true"  estimates  are  required, 
division  by  this  value  is  necessary. 

Another  high-pass  filter  used  in  the  analysis  was  the  difference 
filter  y(t)  =  x(t)  -  0.6x(t  -  l)  which  does  not  atte  'xate  the  low  frequencies 
as  much  as  the  previous  one.  The  spectral  estimates  are  multiplied  by 
(l,2  -  3.6cos2nf)  as  a  result  of  using  this  filter. 


Spectrum  Representation  and  Significance  Tests 

The  variance  spectra  of  the  horizontal  components  have  been  computed 
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using  normalized  data  and  the  Banned  estimates.  These  are  represented  In  the 
simplest  form:  P(f)  ve.  f.  The  trend  in  geophysics  has  been  to  use  logarithmic 
repres entat ion  with  the  spectral  estimate  Iff)  as  the  ordinate  and  log  f  as 
the  abscissa;  also,  f  4  F(f)  vs.  log  f  is  sometimes  used  so  as  to  represent 
the  variance  oy  the  area  between  the  frequencies. 

In  the  spectrum  representation,  the  abscissa  is  the  frequency, 
labeled  from  0.00  to  0.50  cycles  per  data  Interval.  The  data  Intervals  in  the 
real  data  are  2  seconds,  8  seconds,  and  16  seconds. 

If  a  data  set  consisting  of  N  random  numbers  is  used  (N  large),  the 
spectrum  obtained  is  said  to  be  one  of  white  n< ise,  for  the  total  variance  is 
distributed  non- preferentially  (theoretically,  equally)  among  all  the  lines. 
Therefore,  dividing  the  line  variances  by  the  total  variance  as  follows 

for  k  =*  1,  2,  . ..,  M-l 

(3.20) 

for  k  »  0  and  M 

the  "line  ratios"  are  obtained.  The  value  of  unity,  then,  should  be  the  theo¬ 
retical  value  for  white  r^ise.  Thus  in  the  spectrum  representation,  the  ordi¬ 
nate  values  represent  the  ratio  of  the  spectral  estimate  (line  variance)  and 
the  line  variance  from  a  set  of  random  data. 

Pig.  3.2  shows  a  spectrum  of  a  set  of  500  random  numbers  with  a  mean 
value  of  zero  and  a  sample  variance  of  5.  The  spectrum  in  (e)  shows  some 
values  below  the  theoretical  value  of  1.0  and  some  values  above.  Clearly,  in 
order  for  the  null  hypothesis  of  white  noise  to  be  meaningful,  confidence  limits 
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are  needed. 


The  white  noise  hypothesis  made  in  this  analyses  is  that  the  spectral 
estimates  obtained  from  the  wind  data  and  those  obtained  from  a  set  of  random 
numbers  are  the  seme.  Two  confidence  limits  are  computed)  namely,  the  a  priori 
and  the  a  posteriori.  As  the  names  imply,  the  former  makes  the  hypothesis  that 
the  estimates  or  line  ratios  to  be  computed,  will  have  significantly  higher 
peaks  at  some  expected  frequencies,  whereas  the  latter  name  implies  that  the 
spectrum  has  been  computed  and  the  peaks  will  be  found.  The  a  posteriori 
concept  ie  explained  by  Eddy,  et  al  (1968,  pp.  26I-65). 

The  confidence  limits  for  white  3pectra  were  found  to  be  as  follows  1 


95# 

1.99  a  posteriori 

1.50  a  priori 

mean 

1.00 

5# 

0.60  a  priori 

0.37  a  posteriori 

These  were  obtained  with  a  stochastic  model  which  consisted  of  500  spectrum 
analyses  of  500  sets  of  random  numbers.  The  procedure  employed  is  similar  to 
that  of  a  Monte  Carlo  simulation  model.  The  details  are  given  both  by  Duchon 
(1968,  Chap.  V)  and  Almazan  (1968,  Chap.  VII). 

In  spectral  analysis  of  geophysical  data  such  as  that  of  horizontal 
wind  components,  large  peaks  in  the  low  frequency  range  of  the  spectrum  and 
very  iow  values  at  the  high  frequencies  are  found  whenever  there  is  a  large 
autocorrelation  value  at  lag  one,  that  is,  a  large  degree  of  correlation 
between  successive  values,  (it  will  be  shown  that  this  pattern  is  found  in 
all  cases  except  under  extremely  turbulent  conditions.)  This  particular 
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feature  may  be  analyzed  for  statistical  significance  by  ass  urn!  ■  g  u  red  noise 
hypothesis;  that  is ,  assuming  that  the  data  are  the  same  as  those  obtained  by 
a  first-order  linear  Markov  process. 

Gilman,  Fuglister,  and  Mitchell  (1963)  have  derived  the  following 
equation  for  determining  a  family  of  red  noise  curves; 

Iv  -  - - o  (3.21) 

1  -  2p  hn/M  +  p“ 

where  p  is  the  lag  1  autocorrelation  value,  h  is  the  harmonic  and  M  is  the 
maximum  lag.  Note  that  if  p  =  0,  which  iB  a  zero  correlation  between  successive 
data  points,  the  spectrum  is  that  of  white  noise.  This  is  shown  in  Fig.  3*3, 
which  has  normalized  variance  or  line  ratios  as  the  ordinate  values. 

For  the  red  noise  hypothesis,  the  95  percent  confidence  limit  to  be 
shown  in  the  spectra  of  lateral  and  longitudinal  data  were  found  by  "chi- 
square  over  degrees  of  freedom",  or  X  /v.  The  effective  degrees  of  freedom 
for  the  Hanned  spectra  may  be  found  by 

v  ~  gN-  ^  ■**/-  (3-22) 

where  N  is  the  number  in  the  data  set  ( 5C0)  and  M  the  maximum  lag.  This 
equation  for  finding  effective  degrees  of  freedom,  Eq.  (3.22),  is  explained 
in  Blackman  (1966),  and  Eddy,  et  al  (1963,  Chap.  ljj). 

The  values  of  coherence  computed  by  Eq.  (3. 19)  were  tested  by  use 
of  a  table  prepared  by  Mitchell  (l 966)  by  interpolation  of  values  found  by 
Amos  and  Xoopsmans  (1963).  The  95  percent  significance  level  is  found  to  be 
O.283  for  the  Banned  spectra  of  500  values  used  in  horizontal  wind  component 
data  cross-spectra.  Therefore,  any  value  of  0.283  or  higher  was  oaken  to  be 
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Line  Ratio  or  Normalized 


Fig.  3*3  Spectra  of  Red  Noise  for  Lag  X  Autocorrelation 
Values  of  0.0,  0.5;  and  0.8 


statistically  significant. 

In  nearly  all  of  the  distributions ,  the  skewness  statistic  and 
hurt  os is  statistic  were  small  values.  Testing  the  values  for  normality  may 
be  accomplished  by  computing  the  standard  error  (one  standard  deviation)  for 

l/2 

the  skewness  approximated  by  (6/U)  '  (Brooks  and  Carr others,  1953)  where  N 
is  the  number  of  observations.  The  standard  error  for  the  kurtosis  may  be 
roughly  approximated  by  (2k/N)  '  .  These  assume  that  the  N  observations  are 
Independent,  which  is  not  the  case  in  most  of  the  data  samples.  But  the 
autocorrelation  functions  presents  an  indication  of  the  degrees  of  freedom. 

For  example,  if  the  autocorrelation  value  goes  to  zero  at  lag  50,  then  the 
ratio  N/50  may  he  used  as  the  number  of  independent  observations;  hence,  where 
N-2450,  (^o)li/2=  0.35  and  (1.96  x  0.35)  =  0.686  which  means  that  95  times 
out  of  100  the  skewness  value  shall  be  between  -0.686  and  0.686  in  a  normal 
distribution.  The  tests  are  based  on  the  theory  of  errors  and  only  give  a 
first  guess  as  to  the  normality  of  the  distributions. 

D.  Analysis  of  Data  Samples 

1.  Evening  Data  Samples  Taken  at  the  Be^ch  Tower 
Case  I-HNj  Beach  Nighttime^  Sampling  Period  1000  sec  (2213-2217  CDT) 

The  autocorrelation  functions  for  Case  I-BN  are  shown  in  Fig.  3»^(a) 
and  (b)  for  the  lateral  and  longitudinal  components,  respectively.  The  inte¬ 
gral  scales  for  the  longitudinal  case  may  be  approximated  for  the  1  meter  (m) 
level.  For  the  other  two  levels  and  all  three  levels  for  the  lateral  component, 
the  functions  show  that  perhaps  a  longer  period  is  needed  to  approximate  the 
"scale."  The  functions  for  both  the  lateral  and  longitudinal  components  at 
the  1  and  6.7  m  (henceforth  called  the  "j  mrt  level)  levels  show  rapid  up-and- 
down  motion  or  oscillations  superimposed  on  the  autocorrelation  curve.  This 
may  be  indicative  of  high  frequency  (short  wavelength)  components.  Also,  the 
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function  for  the  longitudinal  component  at  1  m  ii, 
for  500  random  numbers. 


similar  to  the  one  in  Fig. 


3.2 


The  normalized  Haimod  spectra  for  the  lateral  components  are  shown  in 
Fig.  3=5  (a),  (b),  (c),  for  the  1,  7,  and  27  m  levels,  respectively.  Even  at  the 
lowest  levels  most  of  the  variance  is  in  the  low  frequencies,  in  fact,  nearly  60 
precent  of  the  total  variance  is  contained  in  the  zero,  1st,  and  2nd  harmonics* 

Applying  the  null  hypothesis  that  the  peaks  are  those  of  white  noise, 
the  hypothesis  is  rejected  only  for  0.01  cycles/2  sec  at  the  95th  precentile. 

The  only  conclusion  one  may  draw  here  is  that  the  data  are  from  a  process  other 
than  white  noise.  As  the  5  percent  a  posteriori  confidence  limit  (not  drawn), 
most  of  the  high  frequencies  are  rejected  at  both  the  bottom  and  middle  levels. 
The  value  for  the  lag  1  autocorrelation  at  the  1  m  level  is  0.8l;  therefore, 
applying  the  null  hypothesis  that  the  difference  between  the  data  spectrum  and 
a  red  noise  spectrum  is  not  significant,  from  Eq.  (5*2)  and  a  significance 
test*  the  95  percent  curve  is  drawn  as  shown  in  Fig.  3.5. 

It  is  seen  that  at  frequencies  0.00.  0.01,  0.40,  0.4-1,  0.42,  and  0.46 
to  0.50  cycles/2  sec,  the  hypothesis  is  rtjs^i,eu.  This  implies  that  at  certain 
frequencies  the  wind  component  data  may  be  similar  to  that  obtained  by  a  first- 
order  linear  Markov  process,  that  is,  that  the  observation  depends  on  the  pre¬ 
vious  observation. 

From  Taylor’s  hypothesis  discussed  in  Chapter  III  relating  the  time 
t  to  the  length  x,  and  the  mean  wind  U,  t  »  x/D,  (referred  to  by  Munn  (1966) 
as  "frozen  turbulence"),  the  frequencies  may  be  converted  to  wavelengths. 

For  example, 

p  =  0.50  cycles/?  sec  =  0.50  cycles 

1.45  meters/  sec  =  2.90  makers 


(P)’1 


5.8  meters/cycle 


5.8  meters  in  wavelength 
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5  Spectra  of  Lateral  Components  at  j.,  7>  and  27  m 
Levels,  Case  I-EN 
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Applying  the  white  and  red  noise  9*5  percentile  confidence  limits,  the  signi¬ 
ficant  frequencies  at  1  m  are  those  corresponding  to  290  »,  and  7.2  r.  arid  the 
region  from  6.3  m  to  5*0  a  in  length. 


n  the  significant  frequencies j  teat  ie,  '  ic  l&i&C  SpetCc 


between  the  significant  frequencies,  is  not  difficult  to  Isr/  gine  since  the 


basic  hypothesis  in  Chapter  II  is  that  there  is  a  slowly  vai.;tng  component 


with  a  rapidly  varying  frequency  component  superimposed,  Thus ,  the  slowly 
varying  component  is  approximately  290  meters  or  larger  in  wavelength,  and 
the  high  frequency  is  approximately  7  meters  and  shorter,  with  the  longer 
waves  containing  most  of  the  energy. 

At  7  m  the  95  percentile  of  the  red  noise  hypothesis  shows  that  the 
significant  frequencies  are  0,00,  0,01,  0,3?,  C, 33,  0.36,  0.37,  0.45  through 
Q«5Q  cycles [?.  sec.  These,  with  the  exception  of  0.00  cycles/2  see,  correspond 
to  wavelengths  of  3^0  m  and  to  the  12  m  i-oglon  and  to  the  region  between  8.4 
and  7.6  meters. 

At  27  m  for  the  lateral  component,  the  spectrum  shows  only  the  0.00 
and  0.01  (664  m  wavelength)  to  be  significant  at  the  95  percentile  of  the 
white  noise  and  red  noise  hypotheses.  These  two  frequencies  account  for  over 
99  percent  of  the  total  variance  indicating  little,  if  any,  eddy  energy  such 
as  that  found  at  the  two  lower  levels  for  shorter  eddies  or  wavelengths .  This 
may  be  substantiated  by  the  5  percent  level  of  significance  for  white  noise 
(not  plotted)  which  indicates  that  the  values  from  frequencies  0.03  to  G.50 
might  not  be  due  to  white  noise  but  to  some  other  cause. 

Prev  .  '  \ening  the  data  with  a  low-bar.d  pass  filter  given  by  the  differ¬ 
ence  filter  x(t)  =  u(t)  -  0„6u(t  -  l),  the  new  analysis  is  shown  in  Fig.  3*6. 
Again,  the  low  frequency  region  contains  nearly  97  percent  of  the  energy.  A 
difference  filter  of  x(t)  =  u(t)  -  u(t  -  l)  applied  to  the  data  gives  the  auto- 
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Fig*  3*6  Spectrum  Analysis  of  Lateral  Component  at  ?7  m.  Data 
Prewhitened  by  Difference  Filter  u(t )-o.6u( t-1 ) , 

Case  I- HI 

54 


correlation  function  curve  shown  in  Fig.  3»7  and  a  spectrum  indicating  that 
frequencies  0.1*2,  O.Vf,  0.1*8  and  0.50  cycles/2  sec  may  he  significant  at  the 
white  noise  95  percentile  a  posteriori  confidence  hand.  These  correspond  to 
wavelengths  of  15-8,  13*8,  and  13.3  meters.  The  red  noise  hypothesis  is 

not  applied  in  this  case  because  the  gain  function  has  increased  the  middle 
and  high  frequencies  so  much  that  the  noise  level  may  be  amplified.  This  noise 
may  be  attributed  to  the  instruments,  recorders,  digitizers,  and  truncation. 
Duchon  (1968)  discusses  the  problem  of  instrument- induced  red  noise  as  a  func¬ 
tion  of  the  time  lag  and  shows  that  the  longer  the  instrument  time  lag,  the 
redder  the  spectrum.  This  also  depends  on  the  sampling  frequency.  The  con¬ 
clusion  to  be  drawn  from  this  is  that  the  wave  length  of  the  physical  wave  may 
be  much  longer  than  our  sampling  time  will  allow,  i.e.  the  "window"  being  used 
for  the  data  at  that  level  is  too  small.  Therefore,  a  sampling  period  longer 
than  1000  sec  should  be  used. 

The  autocorrelation  functions  of  the  longitudinal  components  showed 
that  these  Eulerian  time  correlations  reached  lower  values  than  the  lateral 
components,  but  only  the  one  for  the  1  m  level  reached  a  value  of  zero. 

The  spectra  of  these  longitudinal  components  are  shown  in  Fig.  3.8 
for  the  three  levels.  The  spectrum  for  the  1  m  level  clearly  shows  the  energy 
distributed  over  the  entire  frequency  range .  The  95  percent  white  noise  a  pos¬ 
teriori  confidence  limit  (1.5)  shows  the  first  3  frequencies  (plus  the  one  at 
0.00)  and  O.36  and  0.37  cycles/2  sec  to  significantly  differ  from  white 
noise.  The  red  noise  95  percent  confidence  limit  further  isolates  the  lower 
frequencies  and  emphasizes  the  importance  of  0. 36  cycles/2  sec  and  adds  to  the 
higher  frequencies  the  l*5th  harmonic.  These  correspond  to  wavelengths  of  290, 
li*5,  and  97  meters  in  the  lower  region  and  in  the  high  region  to  about  8  and  6 
meters.  This  compares  with  the  wavelength  of  290,  7  and  6  meters  for  the 
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Fig.  3,8  Spectra  of  Longitudinal  Components  at  1,  7,  and 
27  m  Levels,  Case  I-EN 


lateral  component  at  the  seme  level.  The  short  wavelength  of  6  m  Is  common  to 
hoth  components,  the  7-8  meter  wavelength  Is  perhaps  common,  but  there  are  more 
wavelengths  present  in  the  longitudinal  component. 

At  the  7  m  level  the  frequencies  which  might  be  significant  for  the 
longitudinal  component  correspond  to  wavelengths  of  380  and  190  meters  in  the 
low  frequency  range  and  to  about  11  to  8  meters  in  the  higher  frequencies. 
Comparing  these  to  the  lateral  values,  the  high  frequency  end  1b  the  same  (8m). 
The  other  lateral  wavelengths  were  380,  12,  and  9-8  meters.  Again,  there  are 
more  wavelengths  present  in  the  longitudinal  case.  If  a  crude  simplified  model 
of  the  physical  phenomenon  were  made  based  on  these  two  levels,  one  could  pic¬ 
ture  an  eddy  or  perturbation  resulting  from  differential  heat  radiation  from 
the  earth,s  surface  with  a  horizontal  dimension  of  200  to  4 00  meters.  As  the 
eddy  traveled  along  the  wind,  mechanical  turbulence  induced  by  the  friction  of 
the  surface  boundary  would  produce  the  shorter  wavelengths  obtained  of  6  to  12 
meters;  furthermore,  more  of  these  "disturbances"  would  prevail  in  the  down¬ 
stream  component.  The  greater  energy  found  at  the  higher  frequencies  in  the 
downstream  case  would  come  from  the  mean  motion,  or  longer  waves. 

At  the  27  m  level  there  is  also  more  energy  at  the  high  frequencies 
for  the  downstream  component.  Prewhitening  the  data  shows  that  the  region 
around  0.40  cycles/2  sec  (l6  meter  wavelength)  might  be  significant,  as  shown 
in  Pig.  3.9- 

The  cross- correlations  of  the  lateral  components  and  of  the  longitu¬ 
dinal  components  are  shown  in  Pig.  3*10*  It  appears  that  the  sampling  interval 
is  too  short,  especially  for  the  lateral  components.  However,  the  lateral 
cross-correlations  suggest  that  the  correlation  between  levels  is  nearly  the 
3ame  for  the  7  to  27  m  and  1  to  27  m  intervals  (about  O.90)  and  the  lowest 
correlation  is  between  1  to  27  m  (0.86).  A  possible  explanation  for  thiB  might 
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Pig.  3.10  Cross-correlation  Functions  of  Components*  Case  I-HJ 
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be  made  with  the  simplified  model  of  the  eddy  just  discussed.  As  this  eddy 
moves  downwind,  the  surface  retards  the  eddy  at  the  lower  layer.  Thus  the 
correlation  between  the  lowest  level  and  the  next  higher  one  should  be  some¬ 
what  less  than  the  correlation  between  7  and  27  meters.  Longitudinally,  the 
correlation  between  the  levels  is  not  clear,  although  it  appears  to  be  small. 

The  correlations  obtained  with  a  longer  sampling  period  are  discussed  in  the 
next  section  for  sample  Case  II-BN. 

The  coherence  function  given  by  Eq.  (3-19)  was  computed  for  the  inter- 
and  intra-level  components.  Fig.  3-H  showB  the  coherence  between  the  u  and  v 
components  at  the  1  m  level.  The  95  percent  a  priori  confidence  limit  is  placed 
at  O.283.  The  highest  coherence  is  shown  to  be  at  the  higher  frequencies, 
suggesting  again  a  semblance  of  "horizontal  homogeneity"  at  the  short  wave¬ 
lengths  which  may  be  attributed  to  mechanical  turbulence. 

The  same  results  are  shown  in  Table  III. 3  for  the  values  equal  to  or 
greater  than  O.283.  The  table  shows  the  highest  coherence  values  at  0.37  and 
0,48  cycles/2  sec  at  1  m,  end  at  0.39  cycles/2  sec  at  the  7  m  level. 


Table  III. 3.  Significant  Values  of 
Coherence,  Case  I- BN 


Frequency 

(cyc/2sec) 

Lateral  Component 
l-7m  7-27m  1 

-27m 

Long.  Component 
l-7m  7- 27m  l-27m 

.00 

.938 

.960 

.988 

.  328  None 

.01 

.925 

•  955 

•  979 

>.283 

.02 

.502 

.759 

.293 

.295 

.03 

.849 

.06 

.431 

.09 

.451 

.11 

.392 

•  35 

•  300 

.40 

.307 

•  4l 

.387 

.340 

•  50 
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The  coherence  between  the  lateral  components  show  high  values  at  the 
low  frequencies  for  all  three  level  intervals,  i.e,,  1  to  7,  7  to  37,  and  1  to 
27  meters.  The  coherence  between  1  and  7  a  give  values  which  are  above  the 
level  of  significance  at  the  0. 40  and  0. 4l  cyclea/2  sec  region.  This  corresponds 
to  4.6  m  and  9.5  m  wavelengths  at  1  and  7  meter  levels,  respectively.  This 
could  be  Interpreted  that  as  a  wavelength  of  5  m  passes  the  1  m  level,  a  9-5  c 
wave  occurs  at  the  7  a  level.  This  suggests  a  broadening  of  the  "eddy"  upward, 
as  speculated  in  the  earlier  model.  At  the  other  intervale,  only  the  low  fre¬ 
quencies  (long  wavelengths)  are  correlated. 

The  coherence  between  the  downstream  components  show  significant 
values  at  0.35  cycles/2  sec  between  the  1  and  7  meter  levels  which  correspond 
to  5*4  and  10,9  meters  in  wavelengths.  These  are  somewhat  longer  than  for  the 
lateral  component  (4.6  and  9-5  m)  which  might  be  explained  as  an  elongation  of 
the  eddies  along  the  wind.  Between  the  7  m  and  27  m  levels  the  significant 
wavelengths  are  9.3  and  16- 2  m,  respectively,  which  suggests  that  the  distur¬ 
bance  reaches  at  least  to  27  m  altitude  and  enlarges,  but  is  only  discernible 
in  the  longitudinal  component.  Between  1  and  27  meters  height  there  is  no 
significant  frequency,  which  means  that  an  occurrence  at  the  bottom  level  was 
not  accompanied  by  a  similar  one  consistently  at  the  top  level. 

A  comparison  of  the  integral  scales  of  the  longitudinal  components 
by  taking  one- quarter  of  the  values  at  the  zero  intercept  shows  that  the  values 
obtained  are  approximately  2,  5,  and  9  (seconds)  for  the  1,  7,  and  27  meter 
levels,  respectively.  This  may  be  interpreted  as  an  increase  in  the  scales 
upward  in  the  longitudinal  component.  Thus,  the  "average  size  of  the  eddies" 
may  he  said  to  increase  upward.  (It  should  be  noted,  however,  that  the  low 
end  of  the  spectrum  is  not  precise  because  of  the  small  "window"  used  in  the 
sampling  period.  Therefore,  these  values  are  only  crude  approximations.) 
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Case  II-BNj  Beach  Nightime,  Sampling  Period  4 000  seconds 

This  sampling  period  includes  the  1QQO  seconds  of  Case  I-BN  discussed 
In  the  previous  section;  specifically,  the  first  1000  seconds  are  the  same. 

The  autocorrelation  functions  for  the  lateral  and  longitudinal  com¬ 
ponents  using  a  sampling  time  of  4000  seconds  are  shown  in  Pigs.  3.12  (a)  and 
(b).  (Since  the  sampling  interval  of  one  saxnple/8  seconds  was  obtained  by 
arithmatically  averaging  four  one - s ample -per-2-seeond  values,  this  is  tantamount 
to  using  a  low  pass  filter,  i.e.,  filtering  out  some  of  the  higher  frequencies.) 
The  patterns  suggest  that  the  sampling  period  is  perhaps  adequate,  because  the 
functions  decrease  as  time  increases  and  attain  or  approach  a  value  of  2ero 
correlation.  The  functions  for  both  the  1  a  and  7  a  level  longitudinal  curves 
also  indicate  high  frequency  eddies  superimposed  on  larger  eddies. 

The  spectra  of  the  lateral  components  for  the  1,  7,  and  27  m  levels 
are  shown  in  Pig.  3*13  (a),  (b),  and  (c),  respectively,  with  the  95  percent 
confidence  limits  for  the  white  noise  a  posteriori  hypothesis  and  the  red  noise 
hypothesis  for  lag  1  autocorrelation.  At  the  first  level  the  frequencies 
which  may  be  significant  with  the  corresponding  wavelengths  using  Taylor's 
hypothesis  are  as  follow J 


0.00  cycle/8  sec 

0.01 

133*0  m 

0.29 

46.1 

0.30 

44.5 

0.39 

34-3 

o.4o 

33*^ 

0.45 

29*7 

0.46 

29.0 

0.50 

26.7 

It  is  seen  tha'  this  longer  sampling  period,  while  being  a  better  window  for 
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Spectra  of  Lateral  Components  at  1,  7>  and  27  m 
Levels,  Case  II- BN 
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computing  the  autocorrelations,  sacrifices  the  smaller  wavelengths.  The  spectra 
at  the  T  m  level  shows  significant  peaks  at  O.OO,  0.01,  and  0.02  and  at  0.25 
cyeles/8  sec  and  also  in  the  frequency  region  of  0.42  to  0.4-7  and  0.50  cycles/8 
sec.  This  corresponds  to  wavelengths  of  1672  and  836  meters,  another  of  67 
meters,  and  also  the  region  between  40  meters  to  33  meters. 

The  lateral  component  at  the  27  m  level  showB  peaks  at  frequencies 
corresponding  to  2632  and  1316  meter  wavelengths  and  126  to  120  meters. 

The  downwind  component  (Pig.  3-14  (a))  for  the  1  m  level  shows 
significant  wavelengths  of  1336,  688,  and  445  meters  for  the  long  waves,  and 
about  35  meters  for  the  short  waves.  The  middle  level  shows  significant  wave¬ 
lengths  of  I672,  836,  557  and  62  meters  and  others  between  42  to  38  meters. 

At  the  top  level  the  significant  wavelengths  are  2632  m,  1316  m,  and  the  82  to 
67  meter  region  and  53  meters. 

The  wavelengths  computed  from  the  significant  wavelengths  are  shown 
in  Table  III. 4  for  the  lateral  and  longitudinal  components.  It  is  seen  that 

Table  III.4t  Wavelengths  Computed  from 
Significant  Frequencies,  Case  II- BN 


1  m 

Lat. 

Long. 

7  m 

Lat. 

Long. 

27  m 

Lat. 

1336 

1336 

1672 

1672 

2632 

2632 

668 

668 

836 

836 

1316 

1316 

445 

67 

62 

45 

35 

4o 

126 

82 

34 

42 

120 

67 

29 

36 

38 

53 

27 

33 

for  a  sampling  period  of  4000  seconds  (about  67  minutes)  the  eddy  appears  to 
be  elongated  along  the  wind  at  the  1  and  7  m  levels  because  the  shorter  wave- 
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Fig.  3.14  Spectra  of  Longitudinal  Components  at  1,  7,  and 
27  m  Levels,  Case  II- BN 
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lengths  are  found  with  the  lateral  component.  This  table  does  not  show  the 
relative  energy  at  these  wavelengths.  However,  the  spectra  show  the  energy 
to  be  much  less  at  high  frequencies  for  the  lateral  component.  Therefore, 
mechanical  turbulence  appears  to  be  the  main  cause  of  these  eddies,  and  in  the 
period  of  just  over  one  hour,  an  apparent  elongation  downstream  of  those  eddies 
is  noticed  at  the  lower  two  levels. 

At  the  shorter  sampling  period  of  1000  see  (16.67  min),  the  smallest 
eddies  computed  were  of  the  same  diameter  at  both  components  at  the  lower  levels. 
This  could  still  be  the  case  because  in  the  stable  regime,  the  smaller  eddies 
could  he  circular  in  the  shorter  time  period  and  the  longer  ones  elongated  down¬ 
stream  with  time. 

At  27  meters,  the  spectra  (Figs.  3*13  (c)  and  3*1^  (=))  a™*  Table  III. 4 
show  the  longitudinal  component  with  more  shorter  wavelengths  with  values  of 
82  to  67  meters  and  53  meters,  and  the  lateral  with  values  in  the  126  to  120 
meters  region. 

The  cross-correlations  between  levels  are  shown  in  Fig.  3.15  for  the 
two  components.  The  correlations  for  the  lateral .components  are  higher  than 
those  of  the  longitudinal  components,  indicating  again  longer  wavelengths  pre- 
uominate  in  the  lateral.  Bu4-  more  important  is  that  during  the  period  of  66.67 
minutes  sampling  time,  the  correlation  between  the  1  m  and  7  ®  levels  is  the 
highest,  with  the  7  m  to  27  m  correlation  the  next  highest,  and  the  smallest 
correlation  between  the  1  m  and  27  at  levels.  This  is  true  for  both  components, 
indicating  that  perhaps  the  eddies  extend  to  at  least  7  meters  but  to  a  lesser 
extent  to  the  27  m  level. 

The  correlation  value  of  0.90  for  the  lateral  components  between  the 
1  and  7  m  levels,  compared  with  a  value  of  0,5^  for  the  downstream  component, 
support  the  assumption  of  mechanical  turbulence  because  this  would  cause  more 
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eddy  "diameters"  to  be  found  downstream  than  between  the  two  levelB  in  the 
lateral.  These  somewhat  "constant"  lateral  diameters  would  give  higher  cor¬ 
relations  . 

The  magnitude  spectra  and  coherence  values  for  the  lateral  components 
are  shown  In  Fig.  3«l6.  The  magnitude  spectrum  is  shown  to  illustrate  that 
the  coherence  may  be  a  better  parameter  to  use  In  the  cross-spectral  analysis 
of  wind  components.  In  both  Fig.  3*16  and  Fig.  3. 17,  the  magnitude  values  do 
show  some  peaks  in  line  variances,  but  these  peaks  are  not  as  prominent  as  the 
coherence  values.  As  expected,  the  increases  in  coherence  values  occur  with 
those  of  the  magnitude  spectrum  since  both  parameters  are  computed  using  the 
real  and  complex  parts  of  the  cross-spectra  as  given  in  Eqs.  (3*1^)  and  (3*15). 
The  magnitude  value  is  the  magnitude  of  the  vector  in  a  complex  plane  in  which 
the  co-spectrum  is  the  real  value  and  the  quadrature  spectrum  Is  the  imaginary 
component . 

In  the  high  frequency  region,  only  0.33  cycles/8  sec  is  significant 
for  the  coherence  values  between  1  m  and  7  meters.  This  corresponds  to  wave¬ 
lengths  of  about  10  and  50  meters.  The  highest  values  of  coherence  for  all 
the  cross- spectra  in  Figs.  3*16  and  3*17  on  the  high  frequency  ends  of  the 
spectra  occur  in  the  downstream  components  between  1  m  and  7  m  at  0.35  cycles/ 
8  sec.  The  wavelengths  here  are  38  and  h8  meters  which  compare  to  kO  and  50 
meters  in  the  lateral.  This  may  perhaps  be  interpreted  as  an  "equilibrium" 
diameter  for  an  eddy,  given  the  conditions  In  terrain,  stability,  and  wind 
profile.  Between  1  m  and  27  m,  the  coherence  value  is  also  the  highest  in  the 
high  frequency  region  of  downstream  wind  as  shown  in  Fig.  3*17  (b).  At  this 
frequency,  the  wavelength  at  27  meters  is  75  m  and  at  the  1  meter  level,  it  is 
38  m,  again  illustrating  that  the  eddy  is  perhaps  discemable  at  the  high  level 
in  the  downstream  component. 
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Frtqutncy  (cychs/  8  stc) 

Pig.  3. IT  Magnitude  Spectra  and  Coherence  Values  for 
Longitudinal  Components,  Case  II— HN 


As  noted  In  the  beginning  of  this  section,  the  autocorrelation 


function  crossed  the  aero  intercept  in  all  data  samples  except  the  7  a  lateral 
component.  Therefore ,  the  integral  scale  may  be  assumed  to  exist  for  the  reali¬ 
sations  investigated.  The  values  obtained  using  the  F(o)  value  of  2q.  (l.24) 
give  approximately  equal  values  of  10  sec  for  all  3  longitudinal  components 
and  values  of  about  5,  6,  and  7  sec  for  the  lateral  components  at  1,  7,  and 
27  meters. 

Case  III- BN i  Beach  Nighttime,  Sampling  Period  4000  seconds 

This  case  covers  the  period  of  66.67  minutes  following  the  period 
of  Case  II- BN.  The  atmosphere  has  became  slightly  more  stable. 

The  autocorrelation  functions  are  show  in  Fig.  3*18.  A  comparison 
of  these  curves  to  those  4000  seconds  earlier  shows  that  the  autocorrelation 
curve  for  the  27  m  level  decreases  more  slowly  with  lag,  indicating  perhaps 
more  energy  in  the  low  frequency  end  of  the  spectrum  and  less  energy  in  the 
high  frequencies.  The  atmosphere  has  therefore  apparently  become  more  stable 
at  this  level.  The  wind  speed  of  this  offshore  flow  (land  breeze)  has  increased 
slightly  at  all  three  levels  from  the  values  one  hour  before. 

The  spectra  for  the  1  m  level  are  shown  in  Fig.  3.19.  As  noted  in 
the  previous  case,  the  longitudinal  component  has  more  relative  variance  in 
the  higher  frequencies.  The  highest  significant  frequency  without  prewhitenlng 
is  0.50  cycles/8  sec  for  the  lateral  and  0.40  cycles/8  sec  for  the  longitudinal 
components,  which  correspond  to  30  and  37  meter  wavelengths,  respectively.  The 
eddy  may  be  said  to  be  elongated  with  the  wind  in  this  case  also. 

At  the  7  m  level  the  spectrum  (Fig.  3.20)  of  the  longitudinal  compo¬ 
nent  show  more  energy  at  lower  frequencies  at  the  high  end  of  the  spectrum 
(0.36  to  0.41  cyclee/8  sec)  than  that  of  the  lateral  component  (0.44  to  0.46 
cycles /8  sec),  which  indicates  longer  waves,  in  general,  longitudinally  (42  m 
47  m). 
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Fig.  3.19  Spectra  cf  Components  at  1  m,  Case  III- EH 
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Fig.  3.^0  Spectra  of  Wind  Components  at  ?  m,  Case  III- BN 
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The  spectra  of  the  components  at  the  27  m  level  are  shown  in  Fig.  3- 21. 


Nearly  all  of  the  energy  is  in  the  long  waves.  Prewhitening  the  data  by  the 
filter  x(t)  =  u(t)  -  u(t  -  l),  the  spectra  of  the  same  components  is  shown  in 
Fig.  3*22.  The  longitudinal  component  is  seen  to  contain  significant  peaks  at 
a  higher  low  frequency,  which  corresponds  to  smaller  wavelength,  than  the  lateral 
component  (758  m  vs  1516  m).  For  the  rest  of  the  spectrum,  both  components 
appear  to  be  about  the  same. 

The  cross -correlation  between  levels  for  the  longitudinal  components 
are  shown  in  Fig.  3*23.  As  in  the  previous  hour,  the  correlation  is  greater 
between  1  and  7  meters  and  least  between  1  and  27  meters.  The  longitudinal 
components  correlation  (not  shown)  also  were  similar  to  those  of  the  previous 
hour.  The  correlations  longitudinally,  however,  were  somewhat  higher  than 
before,  indicating  more  stable  conditions. 

The  highest  value  of  coherence  was  again  between  the  longitudinal 
components  at  1  and  7  m  at  0.34  cycles/8  sec.  Hiis  corresponds  to  wavelengths 
of  44  m  and  57  m  as  compared  to  previous  values  of  38  and  48  meters.  The  in¬ 
crease  in  wind  speed  and  stability  probably  accountE  for  these  longer  wave¬ 
lengths  . 

Approximating  the  scales  of  turbulence  at  F(o),  the  values  obtained 
are  compared  to  values  of  the  previous  hour: 
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Case  IV- BN  1  Beach  Nighttime,  Sampling  Period  8000  seconds 

The  sampling  period  of  8000  sec  or  approximately  2  hours  and  13 
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Fig.  3.21  Spectra  of  Wind  Components  at  27  m,  Case  I1I-EK 
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Fig,  3.23  Croa a -correlation  Functions  of  Longltulinal 
Components,  Case  III- BN 
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minutes  includes  the  two  400C  sec  periods  described  above.  As  before,  the  500 
data  points  were  obtained  by  averaging  consecutive  data  points.  This  time  every 
8  samples  (l  sample  per  2  seconds)  were  averaged,  thus  essentially  using  another 
low- pass  filter  which  further  attenuated  the  high  frequencies.  This  practice, 
while  rather  delicate  because  of  the  decimation  which  may  occur  (Blackman,  1966 
Chap.  10)  is  nevertheless  preferable  in  this  case  since  the  other  recourse  would 
be  to  sample  every  8th  data  point,  which  might  lead  to  serious  aliasing  problems. 

The  autocorrelation  functions  are  shown  in  Fig.  3*24.  These  show  that 
there  is  a  trend  in  the  lateral  component  data  and  more  trend  in  the  longitudinal 
component  data.  (The  cause  will  be  discussed  later  in  Sec.  2.)  If  these  were 
the  only  samples  available  for  analysis,  it  would  be  advisable  to  remove  the 
trend;  e.g.,  by  means  of  a  least  squares  routine. 

The  Bpectra  of  the  components  for  the  lower  two  levels  are  shown  in 
Figs.  3.25  and  3.26.  The  27  m  level  data  was  pre whitened  with  a  high- pass  filter 
in  order  to  show  the  high  frequency  portion.  The  spectrum  of  these  data  is 
shown  in  Fig.  3*27.  The  wavelengths  which  correspond  to  the  significant  fre¬ 
quencies  are  shown  in  Table  HI.  5. 


Table  III. 5.  Wavelengths  Computed  From  Significant 
Frequencies  Case  IV-BN 

Wavelengths  (meters) 


Lat. 

In  T 

Long. 

lat. 

^  m  Long. 

lat. 

1  m  Long. 

2784 

2?84 

3536 

3536 

5536 

5536 

1392 

1392 

1768 

1768 

2768 

2768 

928 

928 

1845 

1845 

101 

104 

1384 

103 

80 

95 

98 

1107 

99 

77 

86 
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198 

923 
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66 

74 
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615 

58 

58 
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Fig.  3.26  Spectra  of  Wind  Components  at  7  m  Level,  Case  IV-BN 


Line  Ratio  Line  Rat  in 


Frequency  ( cycles/ !6 sect 


Fig.  3.27  Spectra  of  Components  at  27  m  Level  After 
Pre whitening  of  Data,  Case  IV- BN 
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In  general,  the  elongation  downstream  ie  still  present  at  the  lover  levels  and 
is  more  evident  at  27  m.  This  does  not  contradict  the  earlier  results  for  the 
4000  sec  sampling  periods  because  the  longer  averaging  time  might  have  smoothed 
out  the  smaller  wavelengths,  or  more  likely,  the  longer  sampling  period  does 
not  detect  the  smaller  waves.  During  this  long  period  of  time,  these  long  wave¬ 
lengths,  which  do  not  necessarily  have  to  be  eddy  diameters  per  ae  but  could 
be  waves  caused  by  the  wind  shear  (or  perhaps  internal  gravity  waves),  show  the 
effects  of  stabilization  of  the  atmosphere. 

2.  Evening  Data  Samples  Taken  at  the  Canal  Tower 

The  time  period  of  2213-0027  CUT  is  the  same  as  the  Beach  period. 

The  land  breeze  is  the  dominant  feature  during  the  entire  period.  Light  winds 
prevailed,  especially  at  the  1  meter  level,  at  times  becoming  calm. 

Although  the  wind  directions  were  the  same  at  both  locations  (northerly 
flow),  the  canal  site  was  under  a  more  stable  regime.  The  profiles  of  the  dry 
bulb  temperature  and  wet  bulb  temperature  are  shown  in  Appendix  E  for  the  canal 
and  for  the  beach  (only  the  period  from  0000  to  0100  CET  are  shown  for  the  canal 
due  to  a  malfunction  of  the  dry  bulb  thermistor  prior  to  midnight) .  The  canal 
air  temperature  is  warmer  than  the  beach,  indicating  that  the  cooling  effect 
has  passed  at  the  canal  site  and  has  progressed  to  the  beach.  (The  temperature 
inversion  shown  at  the  beach  is  due  to  radiative  cooling.)  Subsiding  air  is 
present  at  the  canal  site  causing  the  warmer  temperature.  The  lowest  layer  at 
the  canal  is  still  cooler  because  of  radiative  cooling  of  the  ground  at  0000 
hours.  By  01.00  hours,  the  subsiding  air  has  warmed  all  levels  and  a  more  stable 
regime  is  present.  The  subsiding  air,  being  drier  and  warmer,  explains  the 
temperature  profiles  (both  dry  and  vet)  at  0100  CUT  for  the  canal  site. 

Case  I-CNi  Canal  Nighttime,  Sampling  Period  1000  seconds 

This  case  is  at  the  same  time  period  (2213-2230  COT)  as  Case  I- BN 
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for  the  beach  site.  Th*  mean  wind,  however,  Is  lighter  for  this  location  than 
at  the  beach. 

'Hie  autocorrelation  functions  are  shown  in  Pig.  3.28.  As  in  Case  I- 
2H  for  the  beach,  the  curves  Indicate  longer  waves  laterally  than  longitudinally. 
The  longitudinal  curves  for  the  two  bottom  levels  show  a  rather  sinusoidal  con¬ 
figuration;  this  could  be  attributed  to  the  light  winds  with  occasional  periods 
of  calm.  A  longer  sampling  interval  with  averaged  observations  would  perhaps 
be  better  In  this  case.  In  general,  the  curves  are  similar  to  those  of  Case 
I- BN  (Pig.  3.4). 

The  spectruD'  curves  (not  shown)  for  the  three  levels  were  used  to 
obtain  the  wavelengths  shown  in  Table  III.6.  These  wavelengths  are  associated 

Table  IH.6.  Comparison  of  Wavelengths  at  Beach  and  Canal 
Stations  for  Period  2213-2220  CDT 

Beach  I-BN  Canal  I-BN 

"'Is  '  f  m’ '  ’  &T  m  Tm  “  fm  - srs — 


Let 

Long 

Lat 

Long 

Lat 

Long 

Lat 

Long 

lat 

Long 

Lat 

Long 

290 

290 

380 

380 

664 

66k 

162 

162 

332 

332 

470 

470 

145 

190 

332 

332 

166 

166 

235 

97 

12 

11 

14 

17 

4 

6 

no 

8 

13 

15 

3 

3 

10 

83 

13 

15 

7 

7 

8 

8 

7 

65 

9 

9 

6  6  55 

47 

27 

15 


with  frequencies  whose  valueB  were  statistically  significant.  The  wavelengths 
.'btained  from  Case  I-BN  are  also  shown  for  comparison  purposes.  The  computed 
wavelengths  are  a  function  of  the  wind  speeds;  therefore,  a  comparison  of  the 
smallest  or  largest  values  cannot  be  made,  because  the  winds  are  stronger  at 


the  beach.  (Perhaps  a  point  shown  here  is  that  the  size  of  the  smallest  eddies 
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or  wavelengths  are  the  smallest  ones  detected  by  this  method  of  computation 
using  spectral  estimates.  That  smaller  sizes  were  present  is  shown  by  using 
the  slower  winds  at  the  canal  site,,  i.c;.,  3  meters  at  the  canal  vs  6  meterB  at 
the  beach*  Furthermore,  a  sampling  period  of  one  sample  per  second  could  theo¬ 
retically  give  'wavelengths  half  as  long  than  those  obtained.)  The  comparisons 
which  may  be  made  are  between  the  wavelengths  which  could  be  resolved  at  both 
locations  and  also  the  relative  wavelengths  laterally  and  longitudinally. 

From  the  values  shown,  at  the  lowest  level  the  horizontal  "shape"  is 
circular  at  both  locations.  This  circular  configuration  implied  by  the  values 
extends  upward  at  the  beach  to  27  meters,  but  is  not  quite  as  evident  at  the 
7  m  level  at  the  canal.. 

As  in  Case  I- BN,  a  longer  sampling  period  is  needed  to  obtain  more 
general  characteristics  to  make  a  comparison  of  the  two  locations. 

Cases  II-CN  and  III-CN:  Canal  Nighttime,  Sampling  Periods  1000  seconds 

The  autocorrelation  functions  for  Case  II-CN  is  shown  in  Fig.  3*29- 
Compared  to  those  from  the  beach  location  at  the  same  period  of  2213-2320  COT 
(Case  II-BN)  shown  in  Fig.  3*12,  these  curves  show  more  fluctuations  in  both 
components  and  presumably  more  energy  at  higher  frequencies  longitudinally. 

In  fact.  Case  III-CN  shown  in  Fig.  3.3O  for  a  period  between  2320-0027  COT  is 
more  similar  to  the  one  at  the  beach  in  Case  II-BN  (Fig.  3*12) .  Tills  illustrates 
that  the  land  breeze  regime  which  was  at  the  canal  location  between  the  period 
2213-2320  COT  is  similar  to  the  one  present  at  the  beach  one  hour  later,  or 
from  2320-0027  COT,  as  presupposed  by  the  temperature  profiles.  At  the  canal, 
however,  the  lateral  components  indicated  more  energy  In  the  lower  frequencies, 
perhaps  caused  in  part  by  lighter  winds,  implying  a  more  stable  regime. 

The  wavelengths  computed  from  the  significant  frequencies  of  the 
spectra  (not  shown)  are  given  in  Table  III. 7  for  Case  II-CN  and  Case  III-CN. 
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Table  III. 7.  Wavelengths  Computed  From  Significant 


Frequencies,  Cases  II- CN,  III-CN,  and  IV-CN 

. __  Case  IX- CN.  (2213-2320)-  4000  sec  Wavelengths  (m) _ _ 

1  m 

.  Lat  7  m 

Lat . 27  m 

720 

1324 

1324 

1872 

1872 

360 

360 

672 

672 

936 

936 

240 

240 

448 

448 

624 

624 

33 

180 

53 

331 

103 

468 

20 

19 

51 

265 

53 

375 

32 

58 

75 

10 

Ik 

55 

48 

30 

37 

69 

42 

U  ■  0.90  mps 

U  a  1.68  mps 

U  -  2.34  mps 

37 

Case 

III-CN  (2320-0027)  4000 

sec  Wavelengths  (m) 

Lat  1  m 

Long 

lat  7“ 

Long 

Lat 

.  Long 

92S 

mmm 

"  1572“"“ 

'  lSt2 

2544' ■  ' 

25W 

464 

k6k 

836 

836 

1272 

1272 

309 

309 

48 

557 

848 

848 

22 

232 

42 

60 

170 

127 

24 

19 

34 

47 

57 

116 

36 

55 

102 

U  *  I.16  mps 

U  »  2.09  mps 

U  =  3.18  mps 

88 

Case 

IV-CN  (2213-0027)  8000 

sec  Wavelengths  (m) 

T  .  1  m 

Lat 

Long 

Lat  7m 

Long 

lSl£ 

■renai 

2944 

2944 

4272 

4^72 

808 

808 

1472 

1472 

2136 

2136 

539 

539 

981 

981 

152 

1091 

37 

85 

164 

740 

119 

237 

32 

k3 

155 

103 

100 

225 

32 

73 

63 

95 

214 

66 

159 

60 

95 

U  =>  1.01  mps 

U  »  1.84  mps 

U  =  2.67  mps 

93 

93 


There  are  two  other  differences  between  these  two  periods  at  the  canal.  The 
first  is  that  the  mean  wind  was  somewhat  lighter  during  the  earlier  period  (Case 
II- CN),  The  second  difference  is  that  the  short  waves  of  the  lateral  and  longi¬ 
tudinal  component  are  nearly  eq  for  the  7  m  and  2J  m  samples  of  Case  II-CN, 

whereas  Case  III-CN  one  hour  later  indicates  an  elongation  with  the  wind  as  has 
already  been  observed  at  the  beach  site  at  the  same  sampling  time. 

The  similarity  of  the  canal  period  of  Case  III-CN  with  the  later 
period  of  the  beach  has  already  been  noted  for  the  temperature  profiles  and  the 
autocorrelation  functions.  The  spectra  also  indicate  s4milar  elongation  of  the 
wavelengths  with  the  wind  for  these  two  samples  between  the  two  locations. 

As  for  the  earlier  case  for  the  Canal  it  might  be  deduced  that,  using 
this  method  of  analysis,  light  windE  under  stable  conditions  do  not  indicate 
an  elongation  of  the  eddies  (or  wavelengths).  Instead,  a  circular  pattern  is 
denoted  much  the  same  as  indicated  using  a  shorter  sampling  period  of  1000 
seconds. 

Case  IV-CN:  Canal  Nighttime,  Sampling  Period  8000  seconds 

The  wavelengths  derived  from  the  8000  second  sampling  period,  Case 
IV-CN,  shown  in  TableIH.7  also  indicate  a  similar  circular  pattern;  but  the 
mean  winds  are  light,  while  the  sampling  period  is  long.  Therefore,  at  least 
for  these  samples,  there  appears  to  be  a  connection  between  the  magnitude  of 
the  wind  and  the  indication  of  elongation  with  the  wind  of  the  eddies,  depend¬ 
ing  on  the  sampling  period.  It  should  be  noted,  however,  that  these  cases  of 
light  winds  involve  some  periods  of  calm,  and  are  not  steady  winds  of  small 
magnitude.  These  periods  of  calm  might  bo  the  cause  of  these  differences. 

The  autocorrelation  functions  of  the  components  for  8000  seconds 
sampling  period  are  shown  in  Pig.  3.31.  All  of  the  curves  indicate  a  trend 
in  the  data  samples.  The  trend,  indicated  is  longer  and  more  evident  than 
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Pig.  3«31  Autocorrelation  Functions  of  Wind  Components 

Case  IV- CR 


that  at  the  beach  at  the  same  period  (see  Pig.  3.2^).  The  increasing  wind 
speeds,  more  pronounced  at  the  canal,  might  account  for  most  of  the  trend] 
however,  some  of  the  trend  might  he  due  to  the  gradual  change  in  direction  of 
the  wind. 

3.  Daytime  Data  Samples  Taken  at  the  Beach  Tower 
Case  I-BD!  Beach  Daytime,  Sampling  Period  1000  seconds 

This  sampling  period  covered  the  period  of  1000  CDT  to  1017  COT 
when  the  wind  direction  was  from  the  east.  The  regime  could  not  be  Bald  to  be 
a  sea  breeze  but  "becoming"  a  sea  breeze,  although  the  past  history  of  an  air 
parcel  would  be  maritime  air  modified  as  it  went  across  about  two  hundred  meters 
of  land. 

Autocorrelation  curvea  of  the  components  shown  in  Fig.  3.32  shew  very 
similar  patterns  between  the  lateral  and  longitudinal  components,  suggesting 
perhaps  an  approach  to  horizontal  isotropy  (the  curves  would  have  to  be  exactly 
alike  at  a  given  level  for  isotropic  turbulence).  The  curves  show  an  indication 
of  near  equal  distribution  of  energy  laterally  and  longitudinally  at  a  given 
level;  in  the  stable  case  the  lateral  component  had  most  of  the  energy  in  the 
long  waves  and  the  longitudinal  component  had  more  energy  distributed  over  the 
spectrum. 

The  spectra  of  the  data  for  the  1  m  level  is  shown  in  Pig.  3 . 33 . 

These  are  similar  in  the  low  frequency  region  which  1b  associated  with  wave¬ 
lengths  from  360  meters  to  90  meters  and  the  high  frequency  region  of  wave¬ 
lengths  of  8  meters.  There  are  still  longer  scales  in  the  lateral  component; 
however,  an  approximation  of  the  integral  scale  of  turbulence  by  the  energy  at 
the  zero  frequency  gives  a  value  of  k.6  sec  for  the  lateral  and  a  value  of  3-3 
sec  for  the  downstream  component.  (The  same  conclusion  of  relative  scales 
could  perhaps  have  been  reached  by  comparing  the  two  autocorrelation  curves 
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for  the  1  m  level.)  "Die  local  isotropic  turbulence  la  Implied  because  the 
energy  at  the  high  frequency  end  is  nearly  the  same  for  both  components. 

At  7  m  (Fig*  3* 3^)  above  the  ground,  where  the  mean  wind  is  essentially 
the  same  (l.80  mps  at  1  m  and  1.82  mps  at  7  m),  the  low  frequency  portion  of 
the  spectrum  shows  wavelengths  of  36k  and  182  meters  for  both  components,  and 
about  Ik  meters  on  the  lateral  and  about  18,  11,  and  7  meters  for  the  longi¬ 
tudinal  (downstream)  components. 

The  long  meandering  wind  In  the  lateral  sense  for  the  stable  case  is 
no  longer  present,  for  the  longer  wavelengths  are  essentially  the  same.  As  a 
result  of  the  sun’s  heating  of  the  earth’s  surface,  bouyancy  or  convective 
forces  are  now  contributing  to  the  turbulence,  as  well  at  mechanical  effects. 

This  Is  shown  In  the  "circular"  pattern  of  the  wavelengths  at  the  1  m  level. 

At  the  other  two  levels,  there  are  more  high  frequency  contributions  longitu¬ 
dinally,  as  shown  in  Figs.  3.3k  and  3-35.  To  summarize,  the  approximate  signi¬ 
ficant  wavelengths  are  shown  in  Table  III. 8.  It  appears  that  even  though  the 

Table  III. 8.  Wavelengths  Computed  From  Significant 
Frequencies  Case  I-BD 


I<at  . 

Id  t 

- . Lons 

.  .Lat 

^  ®  Long 

Lat 

27  m  T 

Lour 

360 

360 

36k 

36k 

kkk 

kkk 

180 

180 

182 

182 

222 

222 

120 

120 

lk 

18 

17 

11 

90 

90 

12 

10 

12 

11 

9 

8 

8 

7 

wind' is  from  the  east,  and  has  a  maritime  history,  the  surface  heating  is  enough 
to  create  thermal  turbulence.  Also,  at  this  sampling  period  the  upper  levels 
show  some  elongation  laterally.  The  increase  in  wavelengths  upward  is  to  a  much 
lesser  degree  than  with  the  stable  case. 
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Frequency  (  cycles /2  sec ) 


Fig-  3.31*  Spectra  of  Wind  Components  at  7  m  Level,  Case  I- ED 
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Frequency  (  cycles / 2 sec) 


Frequency  (cycles  /  2  sec ) 


Fig*  3-35  Spectra  of  Wind  Components  at  27  m  Level,  Case  I-BD 
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The  cross- correlation  curves  are  shown  in  Pig.  3.36.  There  is  high 
correlation  of  the  same  components  between  levels,  with  the  lowest  value  obtained 
by  correlating  the  1  m  and  27  m  levels.  The  correlations  between  levels  1  m  and 
7  m  and  between  7  m  and  27  m  are  nearly  the  name.  Furthermore,  the  correlations 
of  the  lateral  components  are  nearly  equal  to  those  of  the  longitudinal  components 
for  the  same  levels;  in  the  stable  atmosphere,  the  longitudinal  correlations  be¬ 
tween  levels  were  lower  that  lateral  correlations.  This  indicates  that  the  con¬ 
vective  type  turbulence  equally  influences  the  horizontal  components,  at  least 
up  to  the  size  of  wavelengths  Implied  by  this  type  of  analysis. 

An  important  feature  shown  in  these  correlations  is  that  the  winds, 
although  veering  with  time,  do  not  veer  simultaneously.  It  is  rather  difficult 
to  determine  from  the  correlation  curves  which  level  is  "ahead"  of  the  other. 

An  examination  of  the  plots  of  the  components  for  this  time  period  does  not 
give  a  clear  picture  of  the  veering  either.  It  was  evident  from  studying  the 
"strip  charts",  such  as  the  one  shown  in  Appendix  A,  that  a  decrease  in  wind 
speed  took  place  during  the  transition  period  from  the  land-breeze  to  the  sea- 
breeze  regfme.  This  occurred  first  at  the  top  level,  then  the  middle  level, 
and  the  bottom  level  was  the  last  to  turn.  The  bivariate  statistics  computed 
for  the  actual  components  (Appendix  F)  show  that  Case  I-BD  from  1000-1017  CDT 
and  Case  lII-BD  from  1005-1022  CDT  have  the  lowest  mean  winds.  These  two  cases 
occurred  during  this  transition  period. 

Case  III- HD!  Beach  Daytime- Sampling  Period  1000  seconds 

These  are  for  the  period  1005-1022  CUT  for  winds  from  an  ESE  direction, 
thus  traversing  si  ..ghtly  less  land  area  than  Case  I-BD. 

The  correlation  functions  are  shown  in  Pig,  3« 37 J  to  a  small  degree, 
longer  waves  laterally  are  indicated  for  the  7  m  level  and  27  m  level.  The 
spectra  for  the  three  levels  ale  given  in  Pigs.  3*38,  3«39>  and  3»^0.  The 
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Frequency  (  cycles /  2  sec) 


Spectra  of  Wind  Components  at  1m  Level,  Case  III-BD 


Line  Ratio  Line  Ratio 


(a  >  Lateral 


32.9 


Frequency  <  cycles/ 2  sec) 


Fig.  3 . 4-0  Spectra  of  Wind  Components  at  27  m  Level,  Case  IH-BD 
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figures  in  Table  III. 9  are  obtained  from  these  spectra. 

Table  III. 9*  Wavelengths  Computed  From  Significant 
Frequencies,  Case  III-BD 


Lat 

1  m  T 

Lena 

lat 

7  m  Lon. 

Lat 

87  ”  Icu 

380 

380 

350 

350 

426 

42 6 

190 

190 

175 

175 

231 

231 

127 

127 

117 

11 

39 

13 

95 

15 

18 

8 

18 

10 

16 

8 

14 

7 

13 

9 

8 

The  bottom  level  shows  a  somewhat  circular  configuration  (local 
horizontal  isotropy)  and  the  other  two  levels  show  a  lateral  elongation  and 
smaller  wavelengths  as  in  Case  I-BD.  The  correlations  also  showed  a  similar 
pattern  of  high  correlations  and  indicate  the  winds  turned  at  different  periods 
for  the  three  levels. 

Case  II- BD:  Beach  Day,  Sampling  Interval  1000  seconds 

In  Case  II-BD,  which  occurred  during  the  period  immediately  following 
Case  I-BD,  the  wind  started  to  increase  in  speed.  The  circular  pattern  in  the 
eddies  extends  to  at  least  the  7  m  level.  The  auto-correlations  and  spectra 
are  shown  in  Figs.  3.41  to  3. hi.  Note  that  even  though  the  auto-correlation 
functions  indicate  longer  waves  laterally,  the  spectra  still  show  enough  energy 
at  the  middle  and  high  frequencies  to  give  equal  wavelengths  across-the-wlnd 
and  downwind.  This  is  probably  due  to  increased  surface  heating  causing  con¬ 
vective  turbulence  to  penetrate  at  least  up  to  27  m. 

Case  IV-BDr  Beach  Daytime,  Sampling  Interval  10C0  seconds 

Approximately  thirty  minutes  later  at  104-7-1104  CDT,  the  winds  have 
veered  to  the  SE  and  are  now  definitely  oi'  maritime  origin.  Because  of  the 
increase  In  surface  heating,  the  integral  scales  are  smaller  than  those  during 
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Pig.  j.hj  Spectra  of  Wind  Components  at  7  m  Level,  Case;  II-BD 


Line  Ratio  Line  Rath 


Pig.  'i.kk  Spectra  of  Wind  Components  at  ?7  m  Level,  Case  II- BD 


the  transition  period  from  a  land-breeze  regime  to  a  definite  aea-breeza  regime. 
Fig-  3*45  shows  the  auto- correlation  curves  for  this  period.  The  trend  is 
rather  pronounced  during  this  period  and  is  due  to  the  rapid  veering  of  the 
wind  and.  not  to  the  wind  speeds  as  the  winds  did  not  increase  appreciably 
during  these  17  minutes.  At  the  1  a.  level,  the  spectrum  shows  more  distribution 
of  thi  energy  among  the  entire  frequency  range,  although  long  waves  still  pre¬ 
dominate  (Fig.  3.46).  An  approach  to  horizontal  local  isotropy  is  implied  by 
the  similar  values  at  the  high  frequencies  of  the  components.  At  7  m,  the 
spectrum  (not  shown)  depicted  equal  wavelengths,  and  at  27  m,  the  lateral 
spectrum  showed  higher  frequencies  were  significant j  i.e.,  a  longitudinally 
elongation  this  time,  instead  of  lateral  as  in  the  previous  three  cases  when 
the  wind  was  in  the  transition  period. 

Case  V-EDi  Beach  Daytime,  Sampling  Period  1000  seconds 

Still  later  (1103-1120  CDT),  the  veering  with  time  of  the  wind  has 
decreased  as  indicated  by  the  trend  in  the  correlation  curves  of  Fig.  3.47. 

The  energy  is  even  more  uniformly  distributed  among  the  spectra  as  shown  in 
Fig.  3.48  for  the  1  m  level  and  a  circular  shape  of  the  eddy  is  indicated  for 
this  level.  At  the  7  m  and  27  m  levels,  an  elongation  with  the  wind  is  indi¬ 
cated  by  the  spectra  in  Figs.  3.49  and  J.50. 

Cases  VI- BD  and  VII- BDj  Beach  Dry  time.  Sampling  Periods  4000  seconds 

These  two  cases  which  covered  the  periods  1000-1107  CDT  and  1103-1210 
CDT  are  different  because  the  ^irst  period  covers  the  transition  from  land 
breeze  to  sea  breeze;  the  latter  period  cover’s  the  sea  breeze  regime  only. 

In  this  latter  the  significant  high  frequencies  were  the  same  for  both  compo¬ 
nents;  i.e.,  equal  wavelengths,  for  t^e  first  two  levels.  Thus,  the  bouyancy 
effect  is  noticeableat  large  wavelengths;  e.g.,  about  43  meters  *t  the  bottom 
level  and  54  meters  at  the  T  meter  height.  The  autocorrelation  for  Case  VII- BD 
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Pig.  3.H8  Spectra  of  Wind  Components  at  1  m  Level,  Ce.Be  V-BD 
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Fig.  j.'l9  Spectra  of  Wind  Components  at  7  m  Level,  Case  V-HD 
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Fig.  3*50  Spectra  of  Wind  Components  at  27  m  Level,  Cbbs  V-BD 
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and  the  spectra  for  the  two  lower  levels  are  shown  in  Figs.  3-51  to  3.5;> 

The  cross-correlations  between  levels  were  very  high  using  this 
longer  sampling  period,  as  shown  in  Pig.  3,54  for  Case  VI- HD.  This  implies 
that  the  eddies  or  convective  disturbance  extends  to  at  least  27  m  both  down¬ 
wind  and  across  the  wind. 

4.  Bay  time  Data  Samples  taken  at  the  Canal  Tower 

In  the  Investigation  of  the  samples  from  the  beach  site,  it  was  found 
that  the  sea  breeze  phenomena  was  observed  to  cause  a  gradual  veering  of  the 
winds  with  time.  In  this  manner,  the  sea-breeze  may  be  said  to  travel  inland 
much  the  same  way  as  the  land-breeze  traveled  toward  the  sea. 

Because  of  this  gradual  movement  inland,  the  canal  site  observations 
did  not  register  the  same  change  in  atmospheric  regimes  until  a  later  time. 

From  the  data  obtained,  the  winds  at  the  beach  veered  to  a  sea-breeze  regime 
between  approximately  1005  to  1015  COT,  whereas  the  canal  wind  changed  between 
1155  to  1210  COT,  a  difference  of  nearly  two  hours  between  the  two  locations 
4.8  kilometers  (3  miles)  apart.  Therefore,  not  all  of  the  data  samples  selected 
for  comparison  overlap  In  time. 

Case  I-CDi  Canal  Daytime,  Sampling  Period  1000  seconds 

This  case  covered  the  period  between  1103  and  1120  COTj  a  period 
dominated  by  the  land  breeze  regime  with  winds  from  the  ENE. 

The  autocorrelation  curveB  for  the  two  components  are  shown  in  Fig. 
3.55.  These  indicate  that  1000  seconds  is  an  adequate  sampling  period. 

Another  of  the  features  is  the  similarity  between  the  correlations  of  both  com¬ 
ponents  .  This  is  In  spite  of  the  land  breeze  dominance.  The  affect  of  solar 
radiation  causing  convective  turbulence  is  thus  the  dominant  factor  here  at  all 
three  levels.  Also,  the  similarity  between  the  respective  curves  may  be  indi¬ 
cative  of  some  horizontal  homogeneity  caused  by  this  surface  heating. 

1  On 
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Fig-  3-53  Spectra  of  Wind  Components  at  7  m  Level,  Case  VII- BD 
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The  spectra  for  the  three  levels  are  given  in  Jigs.  3-56,  3* 57,  and 
3.58.  The  similarity  between  the  spectra  of  the  components  at  the  t>ame  level 
indicate  a  close  approach  to  horizontal  local  isotropy j  l.e. ,  that  the  distri¬ 
bution  of  the  energy  among  the  spectrum  is  invariant  with  respect  to  the  coordi¬ 
nate  axis.  The  shape  of  the  eddies  for  this  case  is  nearly  circular  with  the 
top  level  Indicating  slightly  more  energy  at  the  high  frequencies  in  the  lateral 
component,  (elongation  with  the  wind). 

A  comparison  of  this  case  with  the  beach  data  showB  that  the  closest 
comparison  is  with  Case  H-HD  for  the  period  1017-103^  CDT  when  the  winds  were 
already  in  the  sea  breeze  regime.  In  this  case,  for  the  beach  site,  the  shapes 
of  the  eddies  implied  by  the  spectra  were  circular  for  the  three  levels.  Before 
that  time,  the  shapes  indicated  in  Case  I-BD  and  Case  III-BD  were  circular  at 
tne  bottom  and  elongated  laterally  at  the  top  levels.  This  suggest  that  the 
effect  of  surface  roughness  and  thermal  convection  on  the  high  frequency  end 
of  the  spectrum  are  similar.  For  example,  Case  I-BD  has  lighter  offshore  winds 
and  stronger  thermal  convection  whereas  Case  II- ED  has  stronger  onshore  vinds 
and  less  intense  thermal  convection. 

The  crosB-correlation  curves  between  the  same  components  at  different 
levels  is  shown  in  Fig.  3«59*  The  correlations  are  nearly  equal  for  both  com¬ 
ponents  as  expected,  and  the  highest  correlation  is  between  the  1  m  level  and 
the  7  a  level.  The  correlations  between  the  1  m  and  27  m  levels  and  the  7  a 
and  27  a  levels  are  nearly  equal.  Thus  the  eddy  caused  by  convection  may  be 
Bald  to  be  more  evident  between  1  m  and  7  a,  than  between  7  a  and  27  meters. 

The  curves  also  show  that  the  veering  of  the  winds  implied  by  the 
autocorrelation  curves  occurs  with  one  level  ahead  of  the  others,  as  was  the 
case  at  the  beach.  But  these  cross-correlations  also  show  that  the  veering  of 
the  windu  at  the  three  levels  was  not  as  organized  as  at  the  beach,  i.e.,  there 
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?ig*  3 • 56  Spectra  of  Wind  Components  at  1  n  Level,  Case  I-CD 
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were  more  fluctuations  at  the  canal  during  the  period  sampled.  Again,  this 
could  be  attributed  to  the  heat  convection. 

Case  II- CDi  Canal  Daytime,  Sampling  Period  1000  seconds 

This  case  occurred  during  1155-1212  COT,  or  at  noon.  During  this 
time  the  transition  from  a  land  breeze  to  a  sea  breeze  occurred.  As  in  the 
case  studied  at  the  beach  site,  the  transition  was  characterized  by  a  decrease 
In  wind'  speeds  at  all  three  levels. 

The  autocorrelation  curves  and  spectra  for  the  components  at  the 
three  levels  are  shown  In  Figs .  3.6O  to  3.63.  More  energy  in  the  longitudinal 
long  waves  is  indicated  by  the  correlations}  the  spectra  lend  support  to  this. 

Hie  cross-correlation  curves  of  the  components  (Fig.  3*64)  show 
higher  correlations  between  the  downstream  components  than  between  the  lateral 
components,  further  supporting  the  indication  of  longer  waves  downstream  at 
this  sampling  interval  and  at  these  low  wind  speed  values. 

Case  III- CD 1  Canal  Daytime,  Sampling  Period  1000  seconds 

About  thirty  minutes  after  the  change  to  a  sea  breeze  flow,  the 
stronger  wind  was  from  a  line  normal  to  the  coast  line;  thus,  a  flow  of  mari¬ 
time  origin  which  had  traversed  4.3  km  of  land  during  mid-day.  The  autocorre¬ 
lation  curves  and  spectra  are  shown  in  Figs.  3.65  bo  3.68.  At  the  1  m  level 
the  spectrum  (Fig.  3.66)  shows  a  "circular"  pattern  at  the  higher  frequencies, 
corresponding  to  about  7  meters  in  wavelength  and  indicated  local  isotropy. 
The  two  higher  levels  show  a  slight  elongation  downstream,  e.g.,  the  smaller 
wavelengths  at  these  levels  were  about  10  m  laterally  and  11  m  longitudinally. 

A  comparison  of  these  with  the  beach  site  (Cases  IY-BD  and  V-ED) 
shows  that  a  similar  slight  downstream  elongation  also  occurred  after  the  sea 
breeze  was  in  effect.  Prior  to  that  a  lateral  pattern  was  found  in  both  loca¬ 
tions,  for  the  7  o  and  27  m  levels  at  the  beach  (Cases  I-BD  and  HI-BD)  and  at 
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Fig.  3.61  Spectra  of  Wind  Components  at  1  m  Level,  Case  II-CD 
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Fig.  3*62  Spectra  of  Wind  Components  at  7  m  Level,  Case  II-CD 
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Pig.  3.63  Spectra  of  Wind  Components  at  27  m  Level,  Case  II-CD 
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3.66  Spectra  of  Wind  Components  at  1  m  Level,  Case  IH-CD 
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Fig.  3*^7  Spervtra  of  Wind  Components  at  7  m  Levels  Case  III- CD 
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Fig.  3-68  Spectra  of  Wind.  Components  at  27  m  Level,  Case  III-CD 
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the  27  m  level  at  the  canal  (Case  III- CD). 

Cases  IV- CD,  V-CD,  and  VI-CDt  Canal  Daytime,  Sampling  Periods  lOOO  seconds 

These  cases  for  1103-1210  COT,  1155-1302  COT,  and  1210-1317  COT  all 
included  sane  of  the  transition  period.  As  in  the  case  of  the  components  of 
the  beach  study,  the  correlation  corves  show  the  trend  for  this  longer  period 
(Figs.  3.69  and  3.70).  Because  of  the  longer  trends  indicated  In  the  data, 
only  general  characteristics  may  be  obtained  from  the  spectra  (Figs.  3.70,  3*71, 
and  3.72).  The  general  features  of  the  1  m  level  spectra  included  a  similar 
distribution  of  energy  among  the  frequencies  for  both  components  at  the  1  m 
level,  which  is  an  indication  of  horizontal  isotropic  turbulence. 

A  comparison  of  these  sampling  periods  with  those  of  the  stable  evening 
cases  indicates  that  these  longer  sampling  periods  do  not  show  a  distinct  pre¬ 
ferred  elongation  with  or  across  the  wind  at  the  two  upper  levels,  such  as  the 
downstream  elongation  of  the  longer  evening  cases.  During  these  longer  sampling 
periods,  however,  the  cross-correlations  between  levels  (Figs.  3*73  tad  3-74) 
were  again  high  both  laterally  and  longitudinally  compared  to  the  evening  cases 
■which  had  high  correlations  only  for  the  lateral  components. 

Case  VII-CDj  Canal  Daytime,  Sampling  Interval  8000  seconds 

This  case  for  the  period  1103-1317  COT  showed  that  the  sampling 
period  was  longer  than  needed  to  do  this  type  of  study.  The  affects  of  the 
veering  and  the  increase  in  wind  speeds  combined  to  give  a  pronounced  trend  in 
the  data  (Fig.  3.75).  This  would  require  a  trend  removal  procedure  before 
further  analysis  could  be  made.  But  Bince  it  has  been  shown  that  a  1000  sec 
sampling  period  was  adequate  for  analysis  of  the  daytime  cases,  the  trend 
removal  was  not  done  in  this  particular  analysis. 
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3*72  Spectra  of  Wind  Components  at  1  m  Level,  Case  VI-CD 
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IV.  SUMMAT  r  AND  CONCLUSIONS 


The  hypothesis  that  there  are  changes  in  the  turbulent  structure  of 
the  wind  caused  by  the  earth’s  surface  characteristics  and  thermal  variations, 
and  that  these  changes  can  be  detected  by  spectral  analysis  has  been  verified 
by  the  wind  observations  taken  at  two  towers  erected  approximately  4.8  km 
apart  near  the  Texas  Gulf  coast.  The  particular  data  samples  included  the 
land-sea-breeze  phenomenon,  which  allowed  samples  to  be  taken  in  onshore  and 
offshore  wind  flow  during  the  course  of  a  fourteen  hour  experiment.  The  wind 
measurements  were  taken  at  three  elevations  at  both  towers.  Previous  studies 
of  the  wind  structure  have  been  done  using  at  most  two  elevations;  moreover, 
these  earlier  studies  have  not  been  made  in  a  land-Bea-breeze  regime. 

It  was  shown  that  by  recording  and  digitizing  the  data  electronically, 
the  spectra  of  the  wind  velocity  components  may  be  used  to  postulate  the  wind 
structure  by  testing  the  variance  spectral  estimates  for  statistical  significance 
with  two  hypothesis;  namely,  white  noise  and  red  noise.  The  spectra  depicted 
in  the  figures  are  not  composite;  but  rather,  each  spectrum  represents  a  data 
sample.  The  confidence  limits  are  shown  in  each  figure.  Furthermore,  partly 
because  of  the  recording  and  digitizing  techniques  employed,  the  data  did  not 
need  prewhitening  prior  to  the  spectral  analysis.  Only  in  very  stable  (thermal) 
cases  was  there  a  need  to  filter  the  data. 

Throughout  the  data  analysis,  the  assumption  has  been  made  that 
Taylor’s  hypothesis  relating  the  mean  wind,  time,  and  space  is  true  for  all 
data  samples  studied  for  1000  sec  periods  and  4000  sec  periods.  No  attempt  was 
made  to  verify  the  validity  of  this  hypothesis-  the  experiment  would  have  re¬ 
quired  many  wind  sensors  placed  close  together  to  obtain  space  correlations 
and  time  correlations.  Panofsky  et  al  (1958)  investigated  the  validity  of 
Taylor’s  hypothesis  using  data  obtained  by  other  experimenters  and  found  a 
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satisfactory  agreement  for  low  levels  (or  intensities)  of  turbulence.  The 
general  requirement  is  that  the  intensity  be  less  than  unity,  as  shown  theo¬ 
retically  by  Gifford  (1956)*  According  to  Lumley  and  Panofsky  (1964),  the 
limited  information  available  indicates  that  Taylor’s  hypothesis  is  generally 
applicable  except  for  low  wave  number  components,  i.e.,  low  frequency  components, 
of  vertical  velocity  fluctuations.  These  were  not  considered  in  this  study. 

The  bivariate  statistics  for  the  components  are  shown  in  Appendix  F 
for  the  daytime  samples,  using  the  actual  winds.  The  same  statistics  for  the 
downstream  and  lateral  components  are  shown  after  the  rotation  of  the  coordi¬ 
nate  axes.  The  information  obtained  from  the  first  is  the  wind  direction 
(approximate)  measured  with  respect  to  a  North-South  coordinate  system,  and 
the  statistical  parameters  of  the  wind  vector  and  wind  components.  Rotating 
the  axes,  of  course,  does  not  change  the  values  of  some  of  the  parameters  such 
as  the  vector  standard  deviation  of  wind  velocity,  the  vector  coefficient  of 
variation,  or  the  wind  standard  deviations  along  the  minor  and  major  axes  of 
the  ellipse,  since  these  are  all  invariant  with  coordinate  rotations.  After 
rotation  of  the  axes,  the  intensity  of  turbulence,  S^/u,  may  be  computed, 
where  is  the  sample  standard  deviation  of  the  downstream  component  and  U 
is  the  sample  mean  wind.  The  intensities  were  all  less  than  unity  as  shown  in 
the  Tables  F.3  and  F.4.  The  valueB  of  the  Beach  Bamples  were  consistently 
lower  than  those  of  the  Canal.  This  is  due  both  to  smaller  values  of  wind 
speeds  and  greater  standard  deviations  of  the  downstream  wind.  The  smaller 
wind  speed  values  were  a  result  of  the  land-breeze  effect  still  observed  in 
the  canal  area  whereas  a  sea-breeze  had  begun  at  the  beach.  The  larger  stan¬ 
dard  deviations  of  the  downstream  components  were  probably  due  to  the  convec¬ 
tive  forces  on  the  relatively  lighter  wlndB  at  the  lower  measurement  level, 
and  in  some  cases,  by  the  somewhat  erratic  pattern  of  the  turning  of  the  winds 
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towards  onshore  flow  at  all  levels  earned  Indirectly  by  the  thermal  convection 
In  the  land  area.  As  illustrated,  by  Case  III-CD  at  the  canal,  even  after  the 
winds  were  normal  to  the  coast  and  with  the  subsequent  increase  in  wind  speeds, 
the  standard  deviations  at  all  three  levels  were  still  high-  giving  intenslty- 
of- turbulence  values  of  0.37,  0.39,  and  0.32  for  the  1  meter,  7  meter,  and  27 
meter  levels,  respectively.  Panofsky  et  al  (op.  cit.)  found  that  Taylor's 
hypothesis  was  valid  up  to  intensity  values  of  at  least  0.26.  An  experiment 
such  as  the  one  proposed  to  test  Taylor’s  hypothesis  using  a  network  of  anemo¬ 
meters  over  an  extensive  area  would  be  needed  to  determine  the  upper  limit  of 
the  Intensity  of  turbulence  for  Taylor’s  hypothesis  to  be  valid.  For  some  of 
the  beach  samples,  however,  the  values  at  all  three  levels  were  well  below  the 
0.26  value,  e.g. ,  Cases  U-ED,  IV-BD,  and  V-BD. 

After  the  computation  of  the  various  statistical  parameters  and  the 
spectral  analyses  on  the  horizontal  velocity  components,  some  of  the  findings 
are  summarized  as  follows 

In  the  stable  evening  cases,  where  offshore  flow  or  "land-breeze" 
occurred,  oi  e  of  the  prominent  features  found  It)  the  long  lateral  meandering 
waves.  It  is  difficult  to  say  whether  mechanical  turbulence  is  the  foremost 
cause  of  these  long  waves  that  appear  on  the  low  frequency  end  of  the  spectrum 
even  though  it  is  probably  the  cause  of  the  turbulence  which  appears  at  the 
high  frequency  end  at  the  lower  two  levels.  This  is  due  to  the  fact  that  waves 
appear  to  be  quite  separated  from  the  rest  of  the  spectra,  i.e.,  generally  only 
the  firBt  one  or  two  harmonics  are  significant.  The  adjoining  harmonics  have 
uegligable  energy,  a3  shown,  for  example,  in  Fig.  3*5,  for  a  1000  sec  sampling 
period  and  Figs.  3.13  and  3*19  Tor  ^000  sec  sampling  periods  for  the  beach 
evening  cases.  Similar  spectra  were  obtained  from  the  canal  data.  Perhaps 
differential  radlational  cooling  over  the  flat  terrain  induced  these  waves. 
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The  longitudinal  component  in  the  stable  cases  also  showed  more 


energy  at  the  low  frequencies  of  the  spectra.  In  the  spectra  of  the  lateral 
component  the  energy  was  restricted  to  the  first  harmonic]  while  in  the  longi¬ 
tudinal  spectra,  the  energy  was  distributed  among  more  of  the  lower  harmonies. 

A  comparison  of  the  two  components  to  obtain  some  geometric  interpre¬ 
tation  at  the  high  end  of  the  frequency  spectra  shows  that  the  smallest  wave¬ 
lengths  detected  by  this  method  were  on  the  order  of  6  meters  for  the  stable 
case  at  the  beach,  and  3  meters  for  the  canal  at  the  1  m  level  for  both  compo¬ 
nents  for  the  1000  second  sampling  period.  The  suggested  configuration  occas- 
sionally  extended,  upward  to  the  27  meter  level  with  increasing  diameters  up  to 
about  10  to  15  meterB. 

An  increase  of  the  sampling  period  to  4000  seconds  indicated  that 
the  smallest  eddies  detected  were  in  the  order  of  25  to  30  meters.  This  longer 
sampling  period  shown  that  these  eddies  were  elongated  downstream.  However, 
in  cases  in  which  calm  winds  occurred  during  the  sampling  periods  such  as  those 
at  the  canal  site,  the  elongation  was  not  detected.  Because  of  these  periods 
of  calm  at  the  canal,  the  results  for  the  stable  cases  should  not  be  as  reliable 
as  the  beach  data  in  which  no  calm  intervals  occurred. 

The  cross-correlations  of  the  same  components  between  levels  showed 
a  higher  correlation  laterally  than  downstream  implying  more  significant  wave¬ 
lengths  along  the  wind.  The  correlations  between  the  1  m  and  7  m  levels  were 
higher  than  between  7  m  and  27  m  or  1  m  and  27  m  for  the  4000  second  periods. 
Thus  the  25  meter  to  40  meter  diameter  eddies  were  more  discernable  at  the 
lower  levels. 

The  1000  sec  sampling  periods  were  found  to  be  inadequate  for  the 
stable  case  to  fully  detect  the  longer  lateral  wavelengths.  A  4000  sec  period 
was  found  to  be  better,  in  spite  of  the  averaging  process  which  caused  high 
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frequencies  to  be  filtered.  However,  for  the  unstable  (daytime)  case,  the 
1000  sec  period  was  found  to  be  satisfactory.  A  sampling  rate  faster  than  1 
sample  per  2  seconds  was  found  to  be  desirable  in  seme  of  the  day  cases  where 
some  possible  aliasing  was  Indicated. 

The  effect  of  Instability  on  the  lateral  component  spectra  was  found 
to  be  an  increase  of  energy  at  the  low  frequencies  and  little  change  in  the 
high  frequencies.  Thus,  it  appears  that  more  waves  were  superimposed  on  the 
long  waves  of  the  stable  case.  Ac  increase  in  wind,  however,  might  have  an 
effect  on  the  high  frequencies,  according  to  Lumley  and  Fanofslty  (196-b),  by 
increasing  the  energy  of  the  high  frequencies.  The  evidence  is  rather  limited 
(e.g.,  Gases  IX- BD  and  III-BD)  to  reach  any  conclusions  from  this  particular 
study.  For  the  longitudinal  direction,  thermal  stability  had  some  effect  on 
the  low  frequencies  but  not  as  much  as  in  the  lateral  direction.  Wind  speed 
should  have  affected  the  high  frequency  end  as  indicated  by  Ely  (1958),  but 
as  in  the  lateral  component,  the  only  information  of  a  stronger  wind  occurred 
with  an  increase  in  solar  radlational  heating.  Therefore,  it  is  difficult  to 
determine  if  the  increase  in  wind  was  the  factor  in  these  increases  in  energy 
at  the  high  frequencies. 

The  daytime  samples  for  the  1000  sec  periods  showed  that  the  smallest 
eddies  detected  (9  to  14  meters)  were  elongated  downstream  at  the  three  levels, 
but  as  the  solar  radiation  increased,  the  eddies  at  the  lowest  level  became 
circular  as  horizontal  local  isotropy  was  approached  (e.g.,  Cases  III-CD,  III-BD, 
IV- BD,  V-BD).  The  eddies  did  not  increase  in  size  appreciably  upward  as  in  the 
stable  case. 

The  cross-correlation  between  components  at  different  levels  for  the 
unstable  cases  showed  high  values  for  the  longitudinal  and  lateral  components, 
implying  distinct  nearly  circular  eddies  caused  by  convection.  Ibe  correlations 
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of  the  same  component  between  1  a  and  27  a  was  consistently  slightly  lower  than 
those  between  1  m  -  7  m,  and  7  m  -  27  m,  again  implying  a  acre  distinct  eddy 
at  the  lower  levels. 

Hie  cross-correlations  also  showed  that  the  winds  did  not  veer  simul¬ 
taneously  into  a  sea-breeze  flow.  The  graphical  plots  indicated  that  the  top 
level  turned  first,  and  the  bottom  level  last.  The  time  Interval,  however,  was 
on  the  order  of  a  few  minutes. 

A  comparison  of  the  two  locations  during  the  daytime  with  onshore 
flow  indicated  that  the  spectra  at  the  beach  for  the  1  m  level  showed  propor¬ 
tionally  less  energy  at  the  low  frequencies  (Pig.  3.h8)  than  the  spectra  at  1 
m  for  the  canal  (-Pig.  3.66).  Ibis  could  be  due  to  wind  of  maritime  origin 
abruptly  crossing  a  rougher  and  warmer  surface  and  thus  having  the  energy  dis¬ 
tributed  among  higher  frequencies  at  the  expense  of  the  longer  waves.  Por  the 
canal  caseB,  which  were  at  a  later  time,  the  increase  in  solar  radiations! 
heating  could  account  for  the  increase  in  energy  found  in  the  spectra.  At 
the  7  meter  height,  the  increase  was  more  pronounced  at  the  canal  Bite  (Pigs. 
3.b9  and  3*67).  At  the  highest  level  of  27  m,  the  comparison  showed  longer 
waves  at  the  beach  than  at  the  canal  (Figs.  3.50  and  3.58).  This  could  be 
interpreted  as  a  loss  in  energy  of  the  long  waves  due  to  convection  over  land. 

The  explanations  thus  given  of  the  possible  reasons  for  the  informa¬ 
tion  obtained  from  the  spectra  imply  that  if  the  same  three  levels  were  measured 
over  the  water,  most  of  the  energy  would  be  in  the  long  waves,  with  the  longest 
waves  at  the  highest  level.  This  hypothesis  would  have  to  be  examined  by  con¬ 
ducting  such  an  experiment. 

Cross- spectral  analysis  between  the  two  towers  for  2  l/2  hour-record 
periods  at  night,  showed  the  highest  correlation  between  the  lateral  components 
at  27  m;  the  next  highest  between  the  lateral  components  at  7  m  and  the  lowest 
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correlations  between  the  lateral  components  at  the  1  m  level,  and  among  all  of 
the  longitudinal  components.  This  should  be  expected  since  the  27  m  lateral 
components  showed  long  waves  to  be  the  predominant  feature.  Becauae  of  the 
mechanical  turbulence  longitudinally,  the  correlations  between  the  towers  were 
small.  The  coherence  values  were  low  for  all  the  frequencies  except  for  the 
lowest  frequency  for  all  the  cross-spectra  between  the  two  towers,  indicating 
that  the  two  time  series  were  similar  in  long  waves  only. 

The  general  features  which  appeared  to  be  significant  using  1000  sec 
periods  and  onshore  flow  at  the  beach  and  offshore  flow  at  the  canal  (daytime) 
were  the  low  correlation  values  between  the  two  towers  but  high  coherence  values 
at  the  higher  frequencies.  Because  of  the  different  regimes  occurring  simul¬ 
taneously  near  midday,  this  could  be  Interpreted  as  a  detection  of  two  types 
of  turbulence ;  mostly  mechanical  at  the  beach  and  the  other  caused  by  surface 
heating  under  light  wind  conditions  at  the  canal.  This  high  coherence  whs  not 
detected  at  the  higher  levels. 

The  other  cross-spectra  between  the  two  towers  for  the  1000  sec  period 
showed  negligible  correlation  values  and  no  significant  coherence  values.  This 
implies  that  two  entirely  different  time  series  were  being  compared;  that  is, 
that  the  physical  phenomena  at  the  beach  location  were  different  from  that  of 
the  canal.  The  only  exceptions  were  as  mentioned  above  and  alBo  for  the  samples 
taken  30  minutes  after  the  sea  breeze  had  begun.  In  these,  there  were  high  co¬ 
herence  values  at  the  7  and  2?  m  levels  for  low  frequency  and  high  frequency 
ranges  for  the  downstream  component.  These  two  samples  were  not  simultaneous, 
but  nearly  two  hours  apart;  a  comparison  Is  therefore  Impossible  to  make. 

The  coherence  was  found  to  be  a  very  sensitive  measure  of  agreement 
between  two  series.  (Bogart  et  al  (1962)  also  found  this  to  be  the  case,  using 
other  types  of  data.)  Since  the  data  were  digitized  electronically,  the  problem 
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of  noise  in  the  high  frequency  region  introduced  by  hand  digitizing  was  not 
present  in  the  coherence  computations. 

Because  of  the  numerous  types  of  data  interpretation  which  can  be 
made  of  spectral  estimates  of  wind  components,  (such  as  "power  laws”),  a 
deliberate  attempt  has  been  made  to  present  the  spectra  in  graphical  form. 
Other  investigators  might  find  these  of  further  use. 

In  order  to  make  better  interpretations  of  the  wind  phenomena,  the 
vertical  component  is  needed.  Although  the  measurement  of  this  vertical  com¬ 
ponent  at  three  levels  would  necessitate  different  instrumentation,  it  is  re¬ 
commended  that  any  future  experiments  include  the  three  wind  components.  The 
seme  experimental  procedures,  data  reduction,  and  spectral  analysis  d -scribed 
in  the  text  could  be  employed  in  the  investigation  of  the  three-component  wind 
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APPENDIX  B 

Samples  of  plots  of  the  x  -  and  y  -  components  of  Velocity  Page 

B.l  Case  IV-BN  Beach  lower  1  m  Level  12/2213-13/0027  COT  160 

B.2  Case  IV-BN  Beach  Tower  7  a  Level  12/2213-13/0027  COT  i6l 

B.3  Case  IV-BN  Beach  Tower  27  m  Level  12/221 3-1 3/0027  COT  162 

B.4  Case  IV- CN  Canal  Tower  1  a  Level  12/2213-13/0027  CDT  16 3 

B.5  Case  IV- CN  Canal  Tower  7  a  Level  12/2213-13/0027  CDT  164 

B.6  Case  IV-CN  Canal  Tower  27  n  Level  12/2213-13/0027  COT  165 

B.7  Case  IV-BD  Beach  Tower  7  m  Level  13/1047-1104  CDT  166 

B.3  Case  IV-BD  Beach  Tower  27  a  Level  13/1047-1104  COT  167 

B.9  Case  IV- CD  Canal  Tower  1  a  Level  13/1103-1210  CDT  168 

B.10  Case  IV-CD  Canal  Tower  7  a  Level  13/1103-1210  CDT  169 

B.ll  Case  IV-CD  Canal  Tower  27  a  Level  13/1103-1210  COT  170 
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B.l  Sample  of  plotB  of  the  x-  and  y-componente  of  velocity 
Case  IV- BN  Beach  Tower  1  m  Level  12/221 3-13/0027  CDT 


of  the  x-  and  y- components  of  velocity 
rer  27  m  level  12/2213-13/0027  COT 


f  the  x-  and  y-components  of  velocity 
rim  Level  12/2213-13/0027  CHF 


of  the  x-  and  y-components  of  velocity 
■er  7  m  Level  12/2213-13/0027  COT 


y- components  of  velocity 
12/2213-13/0027  COT 


;  of  the  x-  and  y-components  of  velocity 
Tower  7  m  Level  13/1047-1104-  COT 


'  plots  of  the  x-  and  y-components  of  velocity 
Beach  Tower  27  m  Level  13/l047_H04  COT 


Fig.  B.g  Sample  of  plots  of  the  x-  and  y- components  of  velocity 
Case  IV- CD  Canal  Tower  1  m  Level  13/1103-1210  COT 


,e  of  the  x-  and  y-componente  of  velocity 
Tover  7  m  Level  13/1103-1210  CD? 


Pig.  B.ll  Sample  of  plots  of  the  x-  and  y-components  of  velocity 
■  Case  17- CD  Canal  Tower  2J  m  Level  13/il03-1210  COT 


APPENDIX  C 

Statistics  of  Daytime  Samples 

(a)  mean 

(t>)  a2  (variance) 

(a)  0  standard  deviation 

(d)  Skevness 

(e)  Kurtosis 

(f)  X-  from  chi-Bquare  Goodness -of- Pit 


Frequency  Distributions  for  Beach  Daytime  Cases 
Based  on  X-  and  Y- components 


CASE  I 


X 

T 

m 

Y 

7  a 

X  Y 

27  m 

X  Y 

Time 

13/1000  -  1017 

mean 

-1.800 

.092 

-1.670  .723 

-2.217 

-.121 

s2 

.241 

.803 

.554  .555 

.389 

•  755 

s 

.491 

.896 

.745  .745 

.624 

.869 

Skewness 

-.14? 

.106 

-.503  .168 

.881 

.270 

Kurt  os  is 

.484 

-.797 

.743  -.317 

6.485 

-.724 

Chi-square 

60.5 

76.4 

235.4  62.2 

75.3 

98.4 

CASE  II 

1 

21 

7  » 

27  m 

X 

Y 

X  Y 

X 

Y 

Time 

mean 

-2.062 

1.587 

-1.341  2.256 

-2.454 

1.230 

s2 

.378 

.501 

.344  .247 

.265 

•356 

s 

.615 

.708 

.587  .497 

•  515 

.597 

Sk  •.  "ueso 

.171 

-.632 

-.275  -1.565 

-.163 

-.852 

Kurtosis 

-.228 

.940 

-.242  2.601 

-.455 

.565 

Chi-square 

20.5 

62.9 

32.6  273.4 

62.2 

124.0 

CASE  III 

1 

m 

7  m 

27  m 

X 

Y 

X  Y 

X 

Y 

Time 


mean 

-1.828 

.524 

-1.368 

1,089 

-2.095 

UJ 

vn 

CD 

s2 

.256 

1.074 

.274 

.910 

.246 

1.061 

s 

.506 

1.036 

•  523 

•  954 

.496 

1.030 

Skewness 

-.043 

-.109 

-.113 

.040 

2.244 

-.215 

Kurtosis 

.308 

-.716 

3.159 

-.970 

15.191 

-1.063 

Chi-square 

40.2 

59-9 

163.1 

94.1 

45.7 

l4o.5 

' 


t 

Frequency  Distributions  for  Beach  Daytime  Cases 

Based  on 

X-  and  Y-components 

> 

CASE  IV 

1  m 

7  a 

27  m 

X 

Y 

X 

Y 

X 

Y 

Time 

mean 

-1.431 

2.163 

-.610 

2.586 

-1.701 

2.290 

2 

s 

.388 

.230 

.195 

.054 

.106 

.041 

s 

.623 

•  479 

.442 

.232 

.326 

.203 

Skewness 

.335 

-.399 

.351 

-.268 

.124 

-.113 

Kurt oe is 

.269 

.089 

.162 

.o4o 

.061 

.026 

Chi-square 

60.1 

79-2 

21.1 

47.0 

35.1 

26.2 

CASE  V 

1  m 

7  m 

27  m 

X 

Y 

X 

Y 

X 

Y 

Time 

1103  • 

•  1120 

(1000 

sec) 

mean 

-1.194 

2-397 

-.237 

2.670 

-1.530 

2.499 

s2 

.424 

.220 

.133 

.074 

.098 

.039 

s 

.651 

.469 

.364 

-.271 

.313 

.196 

Skewness 

-.232 

-.518 

.295 

-.265 

.274 

.369 

Kurt  os  is 

-.599 

.229 

.027 

-.213 

-.251 

.070 

Chi-6  quare 

108.6 

24.1 

27.4 

28.4 

27.9 

CASE  VI 

1  m 

7  m 

27  m 

X 

Y 

X 

Y 

X 

Y 

Time 

1000  • 

-  1107 

(1000 

sec) 

1  mean 

-1.788 

1.432 

-1.193 

2.031 

-2.135 

1.354 

s2 

.278 

.979 

.446 

.802 

.276 

1.163 

s 

.527 

.989 

.668 

.896 

.525 

1.078 

Skewness 

.113 

-.981 

-.635 

-1.289 

-.256 

-.94 8 

Kurtosis 

-.022 

.171 

.843 

.534 

.304 

-.160 

Chi-square 

21.2 

196.6 

2351.4 

462.5 

63-7 

243.3 

s 
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Frequency  Distributions  for  Beach  Daytime  Cases 
Based  on  X-  and  Y~ components 


CASE  VII 


1  m 

X  Y 

X 

7  a 

Y 

27  m 

X  Y 

Time 

1103  -  1214 

(4ooo 

sec) 

mean 

-.942 

2.540 

.265 

2.693 

-.920 

2.773 

s2 

.297 

.177 

.273 

.110 

.272 

.101 

s 

.545 

.420 

.523 

.332 

.521 

.318 

Skewness 

1 

« 

O 

-.154 

.208 

-.331 

.033 

-.037 

Kurtosis 

-.147 

-.271 

-.347 

•  525 

-.390 

-.465 

Chi-square 

29.7 

23.5 

35.1 

37-5 

35.6 

23-7 
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frequency  Distributions  for  Canal  Day  Cases 
Based  an  X-  and  Y~ component a 


CASE  I 

Level 

Comp 

X 

1  m 

Y 

X 

7  m 

Y 

2?  a 

X  Y 

Time 

1103  - 

■  1120  car 

mean 

-1.278 

-.524 

-2.060 

-.501 

-1.889 

-.699 

s2 

.523 

.566 

1.062 

.890 

.676 

.372 

s 

•  72 

.75 

1.03 

•94 

.82 

.61 

Skewness 

0.421 

.221 

.562 

-.145 

•979 

.123 

Kurtosis 

.488 

-.224 

1.395 

-.725 

2.430 

.250 

Chi-square 

53.1 

34.6 

34.4 

62.7 

39.9 

38.5 

CASE  II 

Level 

1  m 

7  m 

27  m 

Comp 

X 

Y 

X 

Y 

X 

Y 

Time 

1155  -  1212  CDT 

mean 

-.600 

.802 

-.572 

1.565 

-.572 

1.603 

S2 

.561 

1.038 

.777 

1.897 

.429 

.793 

S 

.75 

1.02 

.88 

1.38 

.65 

.890 

Skewness 

.693 

-.277 

-.013 

-.522 

.096 

-.159 

Kurtosis 

.431 

-1.110 

-.106 

-.774 

-.513 

-1.029 

Chi-square 

125.8 

227.5 

4o.o 

282.1 

66.2 

208.0 

CASE  in 

Level 

1  m 

7  « 

27  m 

Comp 

X 

Y 

X 

Y 

X 

Y 

Time 

1236  - 

-  1253  CDT 

mean 

-.510 

1.602 

-.855 

2.335 

-.939 

2.161 

s2 

.667 

.436 

1.161 

1.037 

.781 

.625 

s 

.82 

.66 

1.08 

1.02 

.88 

•79 

Skewness 

.092 

-.137 

-.048 

.080 

-.097 

-.098 

Kurtosis 

-.435 

-.428 

-.800 

-.573 

-.465 

-.307 

Chi-square 

19.8 

34.5 

47.3 

56.8 

52.2 

44.7 
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Frequency  Distributions  for  Canal  Day  Cases 
Based  on  X-  and  Y- components 

CASE  IV 


Level 

Comp 

X 

1  m 

Y 

X 

7  » 

Y 

X 

27  m 

Y 

lime 

1103  - 

1210  CD! 

mean 

-1.279 

-.053 

-1.895 

ro 

v< 

-1.725 

.121 

S2 

.529 

.708 

1.196 

1.415 

.800 

1.110 

S 

.727 

.841 

1.094 

1.190 

.894 

1.050 

Skewness 

.396 

.582 

.322 

.639 

.360 

.822 

Kurtosis 

.506 

.162 

-.178 

.149 

-.297 

.101 

Chi-square 

96.6 

73.0 

27.3 

110.8 

48.9 

124.8 

CASE  V 


Level 

Comp 

X 

1  m 

Y 

X 

7  m 

Y 

X 

27  m 

Y 

Time 

1155  - 

1302  COT 

mean 

-.355 

1.470 

-.335 

2.343 

-.634 

2.129 

s2 

.714 

.667 

1.393 

1.226 

.667 

.591 

s 

.845 

.817 

1.180 

1.107 

.817 

.769 

Skewness 

.298 

-.964 

.287 

-.919 

.303 

-.633 

Kurtosis 

-.399 

.772 

-.451 

•  936 

.000 

.001 

Chi-square  34.7  180.0  26.0  11*8.7  50.1  111.7 

CASE  VI 


Level 

Comp 

X 

1  m 

Y 

X 

7  m 

Y 

X 

27  m 

Y 

Time 

1210  - 

1317  COT 

mean 

-.321 

1.572 

-.283 

2.379 

-.629 

2.123 

s2 

.703 

.458 

1.369 

.883 

.660 

.567 

s 

.838 

.677 

1.110 

•  940 

.822 

.753 

Skewness 

.143 

-.712 

1 — 1 
CO 

H 

-.342 

.220 

-.534 

Kurtosis 

-.345 

.878 

-.432 

-.227 

-.110 

-.125 

Chi-square  45.2 


38.O  34.1  60.4  94.5 
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Frequency  Distributions  for  Canal  Day  Cases 
Based  on  X-  and  Y- components 


CASE  VII 


X 

1  m 

Y 

X 

7  m 

Y 

27  m 

X  Y 

Time 

1103  -  1317  COT 

mean 

-.800 

*759 

-1.089 

1.317 

-1.177 

1.122 

s2 

.801 

1.208 

1.857 

2.213 

.988 

1.807 

s 

.895 

1*099 

1.363 

1.488 

.994 

1.344 

Skewness 

.302 

-.191 

.198 

-.125 

.046 

-.177 

Kurtosis 

-.274 

-1.096 

-.434 

-1.0?4 

-.4oo 

-1.268 

Chi-square 

41.1 

125.4 

33.8 

93.0 

47.9 

169.1 
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APPENDIX  D 

Bivariate  Statistics  for  x  and  y  components 
Table  D.l  Beach  Evening 
Table  D.2  Canal  Evening 
Gain  Function 

Fig.  D.l  Gain  Function  for  First  Difference  Filter  y(t)=<x(t)  -  x(t  -  l) 
Computations : 

1.  Gain  Function 

2.  Spectral  Estimate  Multiplier 
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Table  D.l.  Bivariate  Statistics  for  x-  and  y- components 
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Table  E.l.  Bivariate  Statistics  for  x-  and  y- components 
Beach  Evening 


f. 


Table  D.2.  Bivariate  Statistics  of  x-  and  y-components 
Canal  Evening 


Bivariate  Statistics  of  x-  and  y-eomponents 
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Angle  of  Rotation  of  MaJ.  Axle  17-52 


Computation  of  Gain  Function  for  y(t)  =»  x(t)  -  x(t  -  l)j 
The  response  function  for  x(t)  -  x(t  -  l)  is 


H(f)  =  1  -  e 


~i2rrf 


From  complex  variables : 


sjn  TTf 


2rrf 

e  -  e 

21 


multiplying  by  e"2nf, 


-2nf  ,  1  -  e 

e  sin  TTf  = 


-i2rrf 


21 


or 


0,  -inf  .  -12rrf 

2ie  sin  trf  »  1  -  e 


therefore  H(f)  -  2ie’'"Ln'f  sin  trf 


and  the  gain  for- at  ion, 

Gain  (f)  =  2 | sin  nfj 

2.  Spectral  estimate  Multipliers 

The  spectral  estimate  is  multiplied  by 


1  -  e 


■i2nf 


using  filter  y(t)  =  x(t)  -  x(t  -  l) 


therefore 


(l  -  e"12rrf)(l  -  e-i2nf)  =  2  -  2e"12nf> 


2-2  cos  2nf. 
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Beach  Tower 
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Table  P.2,  Bivariate  Statistics  of  x-  and  y- components 
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Table  F.3»  Bivariate  Statistics  of  Wind  Components  After  Axle 
notation  for  Beach  Daytime  Samples 
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Table  F.3.  Bivariate  Statistics  of  Wind  Components  After  Axis 
Rotation  for  Beach  Daytime  Samples 
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Table  P.4.  Bivariate  Statistics  of  Wind  Components  After 
Rotation  for  Canal  Daytime  Samples 
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Table  F.U.  Bivariate  Statistics  of  Wind  Components  After  Axes 
Rotation  for  Canal  Daytime  Samples 
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Table  F.A.  Bivariate  Statistics  of  Windl  Components  After  Axes 
Rotation  for  Canal  Daytime  Samples 
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Table  F.4.  Bivariate  Statistics  of  Wind  C  mponents  After  Axes 
Rotation  for  Canal  Daytime  Samples 
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